/ Shareware Grab Bag / Shareware Grab Bag.iso / 007 / matrix.arc / MATRIX.S

< prev

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 46 72 6f 6d 20 68 70 6c | 61 62 73 21 68 61 6f 21 |From hpl|abs!hao!|
|00000010| 73 65 69 73 6d 6f 21 72 | 6c 67 76 61 78 21 63 76 |seismo!r|lgvax!cv|
|00000020| 6c 21 75 6d 63 70 2d 63 | 73 21 66 72 65 64 20 57 |l!umcp-c|s!fred W|
|00000030| 65 64 20 44 65 63 20 33 | 31 20 31 36 3a 30 30 3a |ed Dec 3|1 16:00:|
|00000040| 30 30 20 31 39 36 39 0d | 0a 52 65 6c 61 79 2d 56 |00 1969.|.Relay-V|
|00000050| 65 72 73 69 6f 6e 3a 20 | 76 65 72 73 69 6f 6e 20 |ersion: |version |
|00000060| 42 20 32 2e 31 30 2e 31 | 20 20 20 20 20 39 2f 32 |B 2.10.1|     9/2|
|00000070| 37 2f 38 33 3b 20 73 69 | 74 65 20 68 70 6c 61 62 |7/83; si|te hplab|
|00000080| 73 62 2e 55 55 43 50 0d | 0a 50 6f 73 74 69 6e 67 |sb.UUCP.|.Posting|
|00000090| 2d 56 65 72 73 69 6f 6e | 3a 20 76 65 72 73 69 6f |-Version|: versio|
|000000a0| 6e 20 42 20 32 2e 31 30 | 20 35 2f 33 2f 38 33 3b |n B 2.10| 5/3/83;|
|000000b0| 20 73 69 74 65 20 75 6d | 63 70 2d 63 73 2e 55 55 | site um|cp-cs.UU|
|000000c0| 43 50 0d 0a 50 61 74 68 | 3a 20 68 70 6c 61 62 73 |CP..Path|: hplabs|
|000000d0| 62 21 68 70 6c 61 62 73 | 21 68 61 6f 21 73 65 69 |b!hplabs|!hao!sei|
|000000e0| 73 6d 6f 21 72 6c 67 76 | 61 78 21 63 76 6c 21 75 |smo!rlgv|ax!cvl!u|
|000000f0| 6d 63 70 2d 63 73 21 66 | 72 65 64 0d 0a 46 72 6f |mcp-cs!f|red..Fro|
|00000100| 6d 3a 20 66 72 65 64 40 | 75 6d 63 70 2d 63 73 2e |m: fred@|umcp-cs.|
|00000110| 55 55 43 50 0d 0a 4e 65 | 77 73 67 72 6f 75 70 73 |UUCP..Ne|wsgroups|
|00000120| 3a 20 6e 65 74 2e 73 6f | 75 72 63 65 73 0d 0a 53 |: net.so|urces..S|
|00000130| 75 62 6a 65 63 74 3a 20 | 6d 61 74 72 69 78 20 6d |ubject: |matrix m|
|00000140| 61 6e 69 70 75 6c 61 74 | 69 6f 6e 20 73 75 62 72 |anipulat|ion subr|
|00000150| 6f 75 74 69 6e 65 20 6c | 69 62 72 61 72 79 20 69 |outine l|ibrary i|
|00000160| 6e 20 43 0d 0a 4d 65 73 | 73 61 67 65 2d 49 44 3a |n C..Mes|sage-ID:|
|00000170| 20 3c 31 31 38 34 40 75 | 6d 63 70 2d 63 73 2e 55 | <1184@u|mcp-cs.U|
|00000180| 55 43 50 3e 0d 0a 44 61 | 74 65 3a 20 57 65 64 2c |UCP>..Da|te: Wed,|
|00000190| 20 32 37 2d 4a 75 6c 2d | 38 33 20 31 34 3a 33 35 | 27-Jul-|83 14:35|
|000001a0| 3a 32 35 20 50 53 54 0d | 0a 41 72 74 69 63 6c 65 |:25 PST.|.Article|
|000001b0| 2d 49 2e 44 2e 3a 20 75 | 6d 63 70 2d 63 73 2e 31 |-I.D.: u|mcp-cs.1|
|000001c0| 31 38 34 0d 0a 50 6f 73 | 74 65 64 3a 20 57 65 64 |184..Pos|ted: Wed|
|000001d0| 20 4a 75 6c 20 32 37 20 | 31 34 3a 33 35 3a 32 35 | Jul 27 |14:35:25|
|000001e0| 20 31 39 38 33 0d 0a 44 | 61 74 65 2d 52 65 63 65 | 1983..D|ate-Rece|
|000001f0| 69 76 65 64 3a 20 54 68 | 75 2c 20 32 38 2d 4a 75 |ived: Th|u, 28-Ju|
|00000200| 6c 2d 38 33 20 30 30 3a | 32 39 3a 30 31 20 50 44 |l-83 00:|29:01 PD|
|00000210| 54 0d 0a 4f 72 67 61 6e | 69 7a 61 74 69 6f 6e 3a |T..Organ|ization:|
|00000220| 20 55 6e 69 76 2e 20 6f | 66 20 4d 61 72 79 6c 61 | Univ. o|f Maryla|
|00000230| 6e 64 2c 20 43 6f 6d 70 | 75 74 65 72 20 53 63 69 |nd, Comp|uter Sci|
|00000240| 65 6e 63 65 20 44 65 70 | 74 2e 0d 0a 4c 69 6e 65 |ence Dep|t...Line|
|00000250| 73 3a 20 34 37 32 0d 0a | 0d 0a 54 68 65 73 65 20 |s: 472..|..These |
|00000260| 6d 69 67 68 74 20 62 65 | 20 6f 66 20 67 65 6e 65 |might be| of gene|
|00000270| 72 61 6c 20 69 6e 74 65 | 72 65 73 74 2e 20 53 6f |ral inte|rest. So|
|00000280| 72 72 79 20 49 20 64 6f | 6e 27 74 20 68 61 76 65 |rry I do|n't have|
|00000290| 20 61 20 6d 61 6b 65 66 | 69 6c 65 0d 0a 66 6f 72 | a makef|ile..for|
|000002a0| 20 69 74 2c 20 62 75 74 | 20 74 68 65 72 65 20 61 | it, but| there a|
|000002b0| 72 65 20 6e 6f 20 63 6f | 6d 70 6c 65 78 20 64 65 |re no co|mplex de|
|000002c0| 70 65 6e 64 65 6e 63 79 | 20 70 72 6f 62 6c 65 6d |pendency| problem|
|000002d0| 73 2e 0d 0a 0d 0a 3a 20 | 52 75 6e 20 74 68 69 73 |s.....: |Run this|
|000002e0| 20 73 68 65 6c 6c 20 73 | 63 72 69 70 74 20 77 69 | shell s|cript wi|
|000002f0| 74 68 20 22 73 68 22 20 | 6e 6f 74 20 22 63 73 68 |th "sh" |not "csh|
|00000300| 22 0d 0a 50 41 54 48 3d | 3a 2f 62 69 6e 3a 2f 75 |"..PATH=|:/bin:/u|
|00000310| 73 72 2f 62 69 6e 3a 2f | 75 73 72 2f 75 63 62 0d |sr/bin:/|usr/ucb.|
|00000320| 0a 65 78 70 6f 72 74 20 | 50 41 54 48 0d 0a 61 6c |.export |PATH..al|
|00000330| 6c 3d 46 41 4c 53 45 0d | 0a 69 66 20 5b 20 24 31 |l=FALSE.|.if [ $1|
|00000340| 78 20 3d 20 2d 61 78 20 | 5d 3b 20 74 68 65 6e 0d |x = -ax |]; then.|
|00000350| 0a 09 61 6c 6c 3d 54 52 | 55 45 0d 0a 66 69 0d 0a |..all=TR|UE..fi..|
|00000360| 2f 62 69 6e 2f 65 63 68 | 6f 20 27 45 78 74 72 61 |/bin/ech|o 'Extra|
|00000370| 63 74 69 6e 67 20 6d 61 | 74 72 69 78 2e 33 6c 27 |cting ma|trix.3l'|
|00000380| 0d 0a 73 65 64 20 27 73 | 2f 5e 58 2f 2f 27 20 3c |..sed 's|/^X//' <|
|00000390| 3c 27 2f 2f 67 6f 2e 73 | 79 73 69 6e 20 64 64 20 |<'//go.s|ysin dd |
|000003a0| 2a 27 20 3e 6d 61 74 72 | 69 78 2e 33 6c 0d 0a 58 |*' >matr|ix.3l..X|
|000003b0| 2e 54 48 20 4d 41 54 52 | 49 58 20 33 4c 20 31 32 |.TH MATR|IX 3L 12|
|000003c0| 2d 4f 63 74 2d 31 39 38 | 32 0d 0a 58 2e 53 48 20 |-Oct-198|2..X.SH |
|000003d0| 4e 41 4d 45 0d 0a 6d 5f | 63 6f 66 61 63 74 6f 72 |NAME..m_|cofactor|
|000003e0| 2c 20 20 6d 5f 63 6f 70 | 79 2c 20 6d 5f 64 65 74 |,  m_cop|y, m_det|
|000003f0| 65 72 6d 2c 20 6d 5f 69 | 6e 76 65 72 74 2c 20 6d |erm, m_i|nvert, m|
|00000400| 5f 6d 75 6c 74 69 70 6c | 79 2c 20 6d 5f 73 6f 6c |_multipl|y, m_sol|
|00000410| 76 65 2c 0d 0a 6d 5f 74 | 72 61 6e 73 70 6f 73 65 |ve,..m_t|ranspose|
|00000420| 20 5c 2d 20 6d 61 74 72 | 69 78 20 6d 61 6e 69 70 | \- matr|ix manip|
|00000430| 75 6c 61 74 69 6f 6e 20 | 73 75 62 72 6f 75 74 69 |ulation |subrouti|
|00000440| 6e 65 73 0d 0a 58 2e 53 | 48 20 53 59 4e 4f 50 53 |nes..X.S|H SYNOPS|
|00000450| 49 53 0d 0a 58 2e 42 20 | 23 69 6e 63 6c 75 64 65 |IS..X.B |#include|
|00000460| 20 3c 6c 6f 63 61 6c 2f | 6d 61 74 2e 68 3e 0d 0a | <local/|mat.h>..|
|00000470| 58 2e 62 72 0d 0a 58 2e | 42 20 63 63 20 22 2a 2e |X.br..X.|B cc "*.|
|00000480| 63 22 20 2d 6c 6d 61 74 | 72 69 78 0d 0a 58 2e 50 |c" -lmat|rix..X.P|
|00000490| 50 0d 0a 58 2e 42 20 64 | 6f 75 62 6c 65 20 6d 5f |P..X.B d|ouble m_|
|000004a0| 63 6f 66 61 63 74 6f 72 | 28 6d 2c 20 69 2c 20 6a |cofactor|(m, i, j|
|000004b0| 29 3b 0d 0a 58 2e 62 72 | 0d 0a 58 2e 42 20 73 74 |);..X.br|..X.B st|
|000004c0| 72 75 63 74 20 6d 61 74 | 72 69 78 20 2a 6d 3b 0d |ruct mat|rix *m;.|
|000004d0| 0a 58 2e 62 72 0d 0a 58 | 2e 42 20 69 6e 74 20 69 |.X.br..X|.B int i|
|000004e0| 2c 20 6a 3b 0d 0a 58 2e | 50 50 0d 0a 58 2e 42 20 |, j;..X.|PP..X.B |
|000004f0| 73 74 72 75 63 74 20 6d | 61 74 72 69 78 20 2a 6d |struct m|atrix *m|
|00000500| 5f 63 6f 70 79 28 6d 29 | 3b 0d 0a 58 2e 62 72 0d |_copy(m)|;..X.br.|
|00000510| 0a 58 2e 42 20 73 74 72 | 75 63 74 20 6d 61 74 72 |.X.B str|uct matr|
|00000520| 69 78 20 2a 6d 3b 0d 0a | 58 2e 62 72 0d 0a 58 2e |ix *m;..|X.br..X.|
|00000530| 50 50 0d 0a 58 2e 42 20 | 64 6f 75 62 6c 65 20 6d |PP..X.B |double m|
|00000540| 5f 64 65 74 65 72 6d 28 | 6d 29 3b 0d 0a 58 2e 62 |_determ(|m);..X.b|
|00000550| 72 0d 0a 58 2e 42 20 73 | 74 72 75 63 74 20 6d 61 |r..X.B s|truct ma|
|00000560| 74 72 69 78 20 2a 6d 3b | 0d 0a 58 2e 50 50 0d 0a |trix *m;|..X.PP..|
|00000570| 58 2e 42 20 73 74 72 75 | 63 74 20 6d 61 74 72 69 |X.B stru|ct matri|
|00000580| 78 20 2a 6d 5f 69 6e 76 | 65 72 74 28 6d 29 3b 0d |x *m_inv|ert(m);.|
|00000590| 0a 58 2e 62 72 0d 0a 58 | 2e 42 20 73 74 72 75 63 |.X.br..X|.B struc|
|000005a0| 74 20 6d 61 74 72 69 78 | 20 2a 6d 3b 0d 0a 58 2e |t matrix| *m;..X.|
|000005b0| 50 50 0d 0a 58 2e 42 20 | 73 74 72 75 63 74 20 6d |PP..X.B |struct m|
|000005c0| 61 74 72 69 78 20 2a 6d | 5f 6d 75 6c 74 69 70 6c |atrix *m|_multipl|
|000005d0| 79 28 6d 31 2c 20 6d 32 | 29 3b 0d 0a 58 2e 62 72 |y(m1, m2|);..X.br|
|000005e0| 0d 0a 58 2e 42 20 73 74 | 72 75 63 74 20 6d 61 74 |..X.B st|ruct mat|
|000005f0| 72 69 78 20 2a 6d 31 3b | 0d 0a 58 2e 62 72 0d 0a |rix *m1;|..X.br..|
|00000600| 58 2e 42 20 73 74 72 75 | 63 74 20 6d 61 74 72 69 |X.B stru|ct matri|
|00000610| 78 20 2a 6d 32 3b 0d 0a | 58 2e 50 50 0d 0a 58 2e |x *m2;..|X.PP..X.|
|00000620| 42 20 73 74 72 75 63 74 | 20 6d 61 74 72 69 78 20 |B struct| matrix |
|00000630| 2a 6d 5f 73 6f 6c 76 65 | 28 61 2c 20 62 29 3b 0d |*m_solve|(a, b);.|
|00000640| 0a 58 2e 62 72 0d 0a 58 | 2e 42 20 73 74 72 75 63 |.X.br..X|.B struc|
|00000650| 74 20 6d 61 74 72 69 78 | 20 2a 61 3b 0d 0a 58 2e |t matrix| *a;..X.|
|00000660| 62 72 0d 0a 58 2e 42 20 | 73 74 72 75 63 74 20 6d |br..X.B |struct m|
|00000670| 61 74 72 69 78 20 2a 62 | 3b 0d 0a 58 2e 50 50 0d |atrix *b|;..X.PP.|
|00000680| 0a 58 2e 42 20 73 74 72 | 75 63 74 20 6d 61 74 72 |.X.B str|uct matr|
|00000690| 69 78 20 2a 6d 5f 74 72 | 61 6e 73 70 6f 73 65 28 |ix *m_tr|anspose(|
|000006a0| 6d 29 3b 0d 0a 58 2e 62 | 72 0d 0a 58 2e 42 20 73 |m);..X.b|r..X.B s|
|000006b0| 74 72 75 63 74 20 6d 61 | 74 72 69 78 20 2a 6d 3b |truct ma|trix *m;|
|000006c0| 0d 0a 58 2e 53 48 20 44 | 45 53 43 52 49 50 54 49 |..X.SH D|ESCRIPTI|
|000006d0| 4f 4e 0d 0a 54 68 69 73 | 20 73 65 74 20 6f 66 20 |ON..This| set of |
|000006e0| 73 75 62 72 6f 75 74 69 | 6e 65 73 20 64 6f 65 73 |subrouti|nes does|
|000006f0| 20 73 65 76 65 72 61 6c | 20 75 73 65 66 75 6c 20 | several| useful |
|00000700| 6f 70 65 72 61 74 69 6f | 6e 73 20 6f 6e 20 6d 61 |operatio|ns on ma|
|00000710| 74 72 69 63 65 73 20 77 | 68 69 63 68 20 61 72 65 |trices w|hich are|
|00000720| 0d 0a 73 74 6f 72 65 64 | 20 69 6e 20 61 20 73 74 |..stored| in a st|
|00000730| 61 6e 64 61 72 64 20 66 | 6f 72 6d 61 74 20 28 64 |andard f|ormat (d|
|00000740| 65 66 69 6e 65 64 20 69 | 6e 20 2f 75 73 72 2f 69 |efined i|n /usr/i|
|00000750| 6e 63 6c 75 64 65 2f 6c | 6f 63 61 6c 2f 6d 61 74 |nclude/l|ocal/mat|
|00000760| 72 69 78 2e 68 29 2e 0d | 0a 53 70 61 63 65 20 66 |rix.h)..|.Space f|
|00000770| 6f 72 20 6d 61 74 72 69 | 63 65 73 20 69 73 20 61 |or matri|ces is a|
|00000780| 6c 6c 6f 63 61 74 65 64 | 20 61 6e 64 20 66 72 65 |llocated| and fre|
|00000790| 65 64 20 62 79 20 74 68 | 65 20 73 75 62 72 6f 75 |ed by th|e subrou|
|000007a0| 74 69 6e 65 73 0d 0a 58 | 2e 49 20 6d 61 6c 6c 6f |tines..X|.I mallo|
|000007b0| 63 0d 0a 61 6e 64 0d 0a | 58 2e 49 20 66 72 65 65 |c..and..|X.I free|
|000007c0| 2e 0d 0a 4d 61 74 72 69 | 63 65 73 20 73 68 6f 75 |...Matri|ces shou|
|000007d0| 6c 64 20 62 65 20 69 6e | 69 74 69 61 6c 69 7a 65 |ld be in|itialize|
|000007e0| 64 20 62 79 20 74 68 65 | 20 6d 61 63 72 6f 0d 0a |d by the| macro..|
|000007f0| 58 2e 49 20 6d 5f 63 72 | 65 61 74 65 2c 0d 0a 61 |X.I m_cr|eate,..a|
|00000800| 6e 64 20 63 61 6e 20 62 | 65 20 64 69 73 70 6f 73 |nd can b|e dispos|
|00000810| 65 64 20 6f 66 20 62 79 | 20 70 61 73 73 69 6e 67 |ed of by| passing|
|00000820| 20 61 20 70 6f 69 6e 74 | 65 72 20 74 6f 20 74 68 | a point|er to th|
|00000830| 65 20 6d 61 74 72 69 78 | 20 74 6f 0d 0a 58 2e 49 |e matrix| to..X.I|
|00000840| 20 66 72 65 65 2e 0d 0a | 58 2e 50 50 0d 0a 58 2e | free...|X.PP..X.|
|00000850| 49 20 6d 5f 63 6f 66 61 | 63 74 6f 72 5c 63 0d 0a |I m_cofa|ctor\c..|
|00000860| 28 6d 2c 20 69 2c 20 6a | 29 20 72 65 74 75 72 6e |(m, i, j|) return|
|00000870| 73 20 74 68 65 20 63 6f | 66 61 63 74 6f 72 2c 20 |s the co|factor, |
|00000880| 61 73 20 74 79 70 65 0d | 0a 58 2e 49 20 64 6f 75 |as type.|.X.I dou|
|00000890| 62 6c 65 2c 0d 0a 6f 66 | 20 65 6c 65 6d 65 6e 74 |ble,..of| element|
|000008a0| 20 28 69 2c 20 6a 29 20 | 6f 66 20 74 68 65 20 6d | (i, j) |of the m|
|000008b0| 61 74 72 69 78 20 67 69 | 76 65 6e 20 62 79 20 69 |atrix gi|ven by i|
|000008c0| 74 73 20 66 69 72 73 74 | 20 61 72 67 75 6d 65 6e |ts first| argumen|
|000008d0| 74 2e 0d 0a 54 68 61 74 | 20 69 73 3a 20 74 68 65 |t...That| is: the|
|000008e0| 20 64 65 74 65 72 6d 69 | 6e 61 6e 74 20 6f 66 20 | determi|nant of |
|000008f0| 74 68 65 20 6d 61 74 72 | 69 78 20 6f 62 74 61 69 |the matr|ix obtai|
|00000900| 6e 65 64 20 62 79 20 64 | 65 6c 65 74 69 6e 67 20 |ned by d|eleting |
|00000910| 72 6f 77 0d 0a 58 2e 49 | 20 69 0d 0a 61 6e 64 20 |row..X.I| i..and |
|00000920| 63 6f 6c 75 6d 6e 0d 0a | 58 2e 49 20 6a 0d 0a 66 |column..|X.I j..f|
|00000930| 72 6f 6d 20 74 68 65 20 | 6d 61 74 72 69 78 0d 0a |rom the |matrix..|
|00000940| 58 2e 49 20 6d 2c 0d 0a | 6d 75 6c 74 69 70 6c 69 |X.I m,..|multipli|
|00000950| 65 64 20 62 79 0d 0a 58 | 2e 62 72 0d 0a 28 2d 31 |ed by..X|.br..(-1|
|00000960| 29 2a 2a 28 69 2b 6a 29 | 2e 0d 0a 58 2e 50 50 0d |)**(i+j)|...X.PP.|
|00000970| 0a 58 2e 49 20 6d 5f 63 | 6f 70 79 5c 63 0d 0a 28 |.X.I m_c|opy\c..(|
|00000980| 6d 29 20 72 65 74 75 72 | 6e 73 20 61 20 64 75 70 |m) retur|ns a dup|
|00000990| 6c 69 63 61 74 65 20 63 | 6f 70 79 20 6f 66 20 74 |licate c|opy of t|
|000009a0| 68 65 20 6d 61 74 72 69 | 78 0d 0a 58 2e 49 20 6d |he matri|x..X.I m|
|000009b0| 2e 0d 0a 58 2e 50 50 0d | 0a 58 2e 49 20 6d 5f 64 |...X.PP.|.X.I m_d|
|000009c0| 65 74 65 72 6d 5c 63 0d | 0a 28 6d 29 20 72 65 74 |eterm\c.|.(m) ret|
|000009d0| 75 72 6e 73 20 74 68 65 | 20 64 65 74 65 72 6d 69 |urns the| determi|
|000009e0| 6e 61 6e 74 20 6f 66 20 | 74 68 65 20 6d 61 74 72 |nant of |the matr|
|000009f0| 69 78 20 6d 20 61 73 20 | 74 79 70 65 0d 0a 58 2e |ix m as |type..X.|
|00000a00| 49 20 64 6f 75 62 6c 65 | 2e 0d 0a 58 2e 50 50 0d |I double|...X.PP.|
|00000a10| 0a 58 2e 49 20 6d 5f 69 | 6e 76 65 72 74 5c 63 0d |.X.I m_i|nvert\c.|
|00000a20| 0a 28 6d 29 20 72 65 74 | 75 72 6e 73 20 61 20 6d |.(m) ret|urns a m|
|00000a30| 61 74 72 69 78 20 28 69 | 66 20 69 74 20 65 78 69 |atrix (i|f it exi|
|00000a40| 73 74 73 29 20 77 68 69 | 63 68 20 69 73 20 74 68 |sts) whi|ch is th|
|00000a50| 65 20 69 6e 76 65 72 73 | 65 20 6f 66 20 74 68 65 |e invers|e of the|
|00000a60| 20 6d 61 74 72 69 78 0d | 0a 58 2e 49 20 6d 2e 0d | matrix.|.X.I m..|
|00000a70| 0a 58 2e 50 50 0d 0a 58 | 2e 49 20 6d 5f 6d 75 6c |.X.PP..X|.I m_mul|
|00000a80| 74 69 70 6c 79 5c 63 0d | 0a 28 6d 31 2c 20 6d 32 |tiply\c.|.(m1, m2|
|00000a90| 29 20 72 65 74 75 72 6e | 73 20 61 20 6d 61 74 72 |) return|s a matr|
|00000aa0| 69 78 20 77 68 69 63 68 | 20 69 73 20 74 68 65 20 |ix which| is the |
|00000ab0| 72 65 73 75 6c 74 20 6f | 66 20 6d 75 6c 74 69 70 |result o|f multip|
|00000ac0| 6c 79 69 6e 67 20 6d 61 | 74 72 69 78 0d 0a 58 2e |lying ma|trix..X.|
|00000ad0| 49 20 6d 31 0d 0a 62 79 | 20 6d 61 74 72 69 78 0d |I m1..by| matrix.|
|00000ae0| 0a 58 2e 49 20 6d 32 0d | 0a 62 79 20 74 68 65 20 |.X.I m2.|.by the |
|00000af0| 63 6f 6e 76 65 6e 74 69 | 6f 6e 61 6c 20 64 65 66 |conventi|onal def|
|00000b00| 69 6e 69 74 69 6f 6e 20 | 6f 66 20 6d 61 74 72 69 |inition |of matri|
|00000b10| 78 20 6d 75 6c 74 69 70 | 6c 69 63 61 74 69 6f 6e |x multip|lication|
|00000b20| 2e 20 54 68 65 20 6e 75 | 6d 62 65 72 20 6f 66 0d |. The nu|mber of.|
|00000b30| 0a 63 6f 6c 75 6d 6e 73 | 20 69 6e 0d 0a 58 2e 49 |.columns| in..X.I|
|00000b40| 20 6d 31 0d 0a 6d 75 73 | 74 20 62 65 20 74 68 65 | m1..mus|t be the|
|00000b50| 20 73 61 6d 65 20 61 73 | 20 74 68 65 20 6e 75 6d | same as| the num|
|00000b60| 62 65 72 20 6f 66 20 72 | 6f 77 73 20 69 6e 0d 0a |ber of r|ows in..|
|00000b70| 58 2e 49 20 6d 32 2e 0d | 0a 58 2e 50 50 0d 0a 58 |X.I m2..|.X.PP..X|
|00000b80| 2e 49 20 6d 5f 73 6f 6c | 76 65 5c 63 0d 0a 28 61 |.I m_sol|ve\c..(a|
|00000b90| 2c 20 62 29 20 72 65 74 | 75 72 6e 73 20 74 68 65 |, b) ret|urns the|
|00000ba0| 20 6d 61 74 72 69 78 20 | 77 68 69 63 68 20 72 65 | matrix |which re|
|00000bb0| 73 75 6c 74 73 20 66 72 | 6f 6d 20 63 6f 6d 70 75 |sults fr|om compu|
|00000bc0| 74 69 6e 67 3a 20 28 61 | 27 5c 20 61 29 2a 2a 28 |ting: (a|'\ a)**(|
|00000bd0| 2d 31 29 5c 20 61 27 5c | 20 62 2e 0d 0a 54 68 69 |-1)\ a'\| b...Thi|
|00000be0| 73 20 69 73 20 75 73 65 | 66 75 6c 20 66 6f 72 20 |s is use|ful for |
|00000bf0| 73 6f 6c 76 69 6e 67 20 | 61 20 73 65 74 20 6f 66 |solving |a set of|
|00000c00| 20 6c 69 6e 65 61 72 20 | 65 71 75 61 74 69 6f 6e | linear |equation|
|00000c10| 73 20 65 78 70 72 65 73 | 73 65 64 20 61 73 20 6d |s expres|sed as m|
|00000c20| 61 74 72 69 63 65 73 3a | 0d 0a 61 5c 20 78 5c 20 |atrices:|..a\ x\ |
|00000c30| 3d 5c 20 62 2e 0d 0a 58 | 2e 50 50 0d 0a 58 2e 49 |=\ b...X|.PP..X.I|
|00000c40| 20 6d 5f 74 72 61 6e 73 | 70 6f 73 65 5c 63 0d 0a | m_trans|pose\c..|
|00000c50| 28 6d 29 20 72 65 74 75 | 72 6e 73 20 74 68 65 20 |(m) retu|rns the |
|00000c60| 6d 61 74 72 69 78 20 77 | 68 69 63 68 20 69 73 20 |matrix w|hich is |
|00000c70| 74 68 65 20 74 72 61 6e | 73 70 6f 73 65 20 6f 66 |the tran|spose of|
|00000c80| 0d 0a 58 2e 49 20 6d 2e | 0d 0a 58 2e 50 50 0d 0a |..X.I m.|..X.PP..|
|00000c90| 58 2e 49 20 6d 5f 63 72 | 65 61 74 65 5c 63 0d 0a |X.I m_cr|eate\c..|
|00000ca0| 28 6d 2c 20 72 2c 20 63 | 29 20 28 77 68 69 63 68 |(m, r, c|) (which|
|00000cb0| 20 69 73 20 61 20 6d 61 | 63 72 6f 2c 20 6e 6f 74 | is a ma|cro, not|
|00000cc0| 20 61 20 66 75 6e 63 74 | 69 6f 6e 29 20 61 6c 6c | a funct|ion) all|
|00000cd0| 6f 63 61 74 65 73 20 73 | 70 61 63 65 20 66 6f 72 |ocates s|pace for|
|00000ce0| 2c 0d 0a 61 6e 64 20 69 | 6e 69 74 69 61 6c 69 7a |,..and i|nitializ|
|00000cf0| 65 73 20 74 68 65 20 6d | 61 74 72 69 78 0d 0a 58 |es the m|atrix..X|
|00000d00| 2e 49 0d 0a 6d 0d 0a 74 | 6f 20 68 61 76 65 0d 0a |.I..m..t|o have..|
|00000d10| 58 2e 49 20 72 0d 0a 72 | 6f 77 73 20 61 6e 64 0d |X.I r..r|ows and.|
|00000d20| 0a 58 2e 49 20 63 0d 0a | 63 6f 6c 75 6d 6e 73 2e |.X.I c..|columns.|
|00000d30| 0d 0a 58 2e 50 50 0d 0a | 58 2e 49 20 6d 5f 76 5c |..X.PP..|X.I m_v\|
|00000d40| 63 0d 0a 28 6d 2c 20 72 | 2c 20 63 29 20 69 73 20 |c..(m, r|, c) is |
|00000d50| 61 20 6d 61 63 72 6f 20 | 75 73 65 64 20 74 6f 20 |a macro |used to |
|00000d60| 61 63 63 65 73 73 20 65 | 6c 65 6d 65 6e 74 0d 0a |access e|lement..|
|00000d70| 58 2e 49 20 28 72 2c 20 | 63 29 0d 0a 6f 66 20 74 |X.I (r, |c)..of t|
|00000d80| 68 65 20 6d 61 74 72 69 | 78 0d 0a 58 2e 49 20 6d |he matri|x..X.I m|
|00000d90| 2e 0d 0a 54 68 65 20 65 | 6c 65 6d 65 6e 74 73 20 |...The e|lements |
|00000da0| 61 72 65 20 61 6c 77 61 | 79 73 20 6f 66 20 74 79 |are alwa|ys of ty|
|00000db0| 70 65 20 64 6f 75 62 6c | 65 2e 0d 0a 58 2e 53 48 |pe doubl|e...X.SH|
|00000dc0| 20 41 55 54 48 4f 52 0d | 0a 46 72 65 64 20 42 6c | AUTHOR.|.Fred Bl|
|00000dd0| 6f 6e 64 65 72 20 3c 66 | 72 65 64 40 75 6d 63 70 |onder <f|red@umcp|
|00000de0| 2d 63 73 3e 0d 0a 58 2e | 53 48 20 22 53 45 45 20 |-cs>..X.|SH "SEE |
|00000df0| 41 4c 53 4f 22 0d 0a 6d | 61 6c 6c 6f 63 28 33 29 |ALSO"..m|alloc(3)|
|00000e00| 2c 20 66 72 65 65 28 33 | 29 2c 20 70 72 69 6e 74 |, free(3|), print|
|00000e10| 66 28 33 29 0d 0a 58 2e | 53 48 20 44 49 41 47 4e |f(3)..X.|SH DIAGN|
|00000e20| 4f 53 54 49 43 53 0d 0a | 4d 6f 73 74 20 72 6f 75 |OSTICS..|Most rou|
|00000e30| 74 69 6e 65 73 20 72 65 | 74 75 72 6e 20 74 68 65 |tines re|turn the|
|00000e40| 20 63 6f 6e 73 74 61 6e | 74 0d 0a 58 2e 49 20 4d | constan|t..X.I M|
|00000e50| 5f 4e 55 4c 4c 0d 0a 74 | 6f 20 72 65 70 6f 72 74 |_NULL..t|o report|
|00000e60| 20 61 6e 20 65 72 72 6f | 72 2e 20 53 6f 6d 65 20 | an erro|r. Some |
|00000e70| 61 6c 73 6f 20 77 72 69 | 74 65 20 6d 65 73 73 61 |also wri|te messa|
|00000e80| 67 65 73 20 74 6f 20 74 | 68 65 20 73 74 61 6e 64 |ges to t|he stand|
|00000e90| 61 72 64 20 6f 75 74 70 | 75 74 20 76 69 61 0d 0a |ard outp|ut via..|
|00000ea0| 58 2e 49 20 70 72 69 6e | 74 66 2e 0d 0a 58 2e 53 |X.I prin|tf...X.S|
|00000eb0| 48 20 46 49 4c 45 53 0d | 0a 58 2f 75 73 72 2f 69 |H FILES.|.X/usr/i|
|00000ec0| 6e 63 6c 75 64 65 2f 6c | 6f 63 61 6c 2f 6d 61 74 |nclude/l|ocal/mat|
|00000ed0| 72 69 78 2e 68 0d 0a 58 | 2e 62 72 0d 0a 58 2f 75 |rix.h..X|.br..X/u|
|00000ee0| 73 72 2f 6c 69 62 2f 6c | 69 62 6d 61 74 72 69 78 |sr/lib/l|ibmatrix|
|00000ef0| 2e 61 0d 0a 58 2e 53 48 | 20 42 55 47 53 0d 0a 45 |.a..X.SH| BUGS..E|
|00000f00| 72 72 6f 72 20 72 65 70 | 6f 72 74 69 6e 67 20 76 |rror rep|orting v|
|00000f10| 69 61 0d 0a 58 2e 49 20 | 70 72 69 6e 74 66 0d 0a |ia..X.I |printf..|
|00000f20| 73 68 6f 75 6c 64 20 62 | 65 20 61 6e 20 6f 70 74 |should b|e an opt|
|00000f30| 69 6f 6e 2e 0d 0a 2f 2f | 67 6f 2e 73 79 73 69 6e |ion...//|go.sysin|
|00000f40| 20 64 64 20 2a 0d 0a 6d | 61 64 65 3d 54 52 55 45 | dd *..m|ade=TRUE|
|00000f50| 0d 0a 69 66 20 5b 20 24 | 6d 61 64 65 20 3d 20 54 |..if [ $|made = T|
|00000f60| 52 55 45 20 5d 3b 20 74 | 68 65 6e 0d 0a 09 2f 62 |RUE ]; t|hen.../b|
|00000f70| 69 6e 2f 63 68 6d 6f 64 | 20 36 34 34 20 6d 61 74 |in/chmod| 644 mat|
|00000f80| 72 69 78 2e 33 6c 0d 0a | 09 2f 62 69 6e 2f 65 63 |rix.3l..|./bin/ec|
|00000f90| 68 6f 20 2d 6e 20 27 09 | 27 3b 20 2f 62 69 6e 2f |ho -n '.|'; /bin/|
|00000fa0| 6c 73 20 2d 6c 64 20 6d | 61 74 72 69 78 2e 33 6c |ls -ld m|atrix.3l|
|00000fb0| 0d 0a 66 69 0d 0a 2f 62 | 69 6e 2f 65 63 68 6f 20 |..fi../b|in/echo |
|00000fc0| 27 45 78 74 72 61 63 74 | 69 6e 67 20 6d 61 74 2e |'Extract|ing mat.|
|00000fd0| 68 27 0d 0a 73 65 64 20 | 27 73 2f 5e 58 2f 2f 27 |h'..sed |'s/^X//'|
|00000fe0| 20 3c 3c 27 2f 2f 67 6f | 2e 73 79 73 69 6e 20 64 | <<'//go|.sysin d|
|00000ff0| 64 20 2a 27 20 3e 6d 61 | 74 2e 68 0d 0a 73 74 72 |d *' >ma|t.h..str|
|00001000| 75 63 74 20 6d 61 74 72 | 69 78 20 7b 0d 0a 09 69 |uct matr|ix {...i|
|00001010| 6e 74 20 6d 5f 72 6f 77 | 73 2c 20 6d 5f 63 6f 6c |nt m_row|s, m_col|
|00001020| 73 3b 0d 0a 09 7d 3b 0d | 0a 0d 0a 73 74 72 75 63 |s;...};.|...struc|
|00001030| 74 20 6d 5f 20 7b 0d 0a | 09 69 6e 74 20 6d 5f 72 |t m_ {..|.int m_r|
|00001040| 6f 77 73 2c 20 6d 5f 63 | 6f 6c 73 3b 0d 0a 09 64 |ows, m_c|ols;...d|
|00001050| 6f 75 62 6c 65 20 6d 5f | 76 61 6c 75 65 5b 31 5d |ouble m_|value[1]|
|00001060| 3b 0d 0a 09 7d 3b 0d 0a | 0d 0a 64 6f 75 62 6c 65 |;...};..|..double|
|00001070| 20 6d 5f 63 6f 66 61 63 | 74 6f 72 28 29 2c 20 6d | m_cofac|tor(), m|
|00001080| 5f 64 65 74 65 72 6d 69 | 6e 61 6e 74 28 29 3b 0d |_determi|nant();.|
|00001090| 0a 73 74 72 75 63 74 20 | 6d 61 74 72 69 78 20 2a |.struct |matrix *|
|000010a0| 6d 5f 63 6f 70 79 28 29 | 2c 20 2a 6d 5f 69 6e 76 |m_copy()|, *m_inv|
|000010b0| 65 72 74 28 29 2c 20 2a | 6d 5f 74 72 61 6e 73 70 |ert(), *|m_transp|
|000010c0| 6f 73 65 28 29 2c 20 2a | 6d 5f 6d 75 6c 74 69 70 |ose(), *|m_multip|
|000010d0| 6c 79 28 29 2c 20 2a 6d | 5f 73 6f 6c 76 65 28 29 |ly(), *m|_solve()|
|000010e0| 3b 0d 0a 0d 0a 23 64 65 | 66 69 6e 65 09 6d 5f 76 |;....#de|fine.m_v|
|000010f0| 28 6d 2c 20 72 2c 20 63 | 29 09 28 28 73 74 72 75 |(m, r, c|).((stru|
|00001100| 63 74 20 6d 5f 20 2a 29 | 6d 29 2d 3e 6d 5f 76 61 |ct m_ *)|m)->m_va|
|00001110| 6c 75 65 5b 72 20 2a 20 | 28 6d 2d 3e 6d 5f 63 6f |lue[r * |(m->m_co|
|00001120| 6c 73 29 20 2b 20 63 5d | 0d 0a 23 64 65 66 69 6e |ls) + c]|..#defin|
|00001130| 65 09 4d 5f 4e 55 4c 4c | 09 28 28 73 74 72 75 63 |e.M_NULL|.((struc|
|00001140| 74 20 6d 61 74 72 69 78 | 20 2a 29 30 29 0d 0a 0d |t matrix| *)0)...|
|00001150| 0a 23 64 65 66 69 6e 65 | 09 6d 5f 63 72 65 61 74 |.#define|.m_creat|
|00001160| 65 28 6d 2c 20 72 2c 20 | 63 29 20 7b 5c 0d 0a 09 |e(m, r, |c) {\...|
|00001170| 09 09 69 66 20 28 28 28 | 69 6e 74 29 28 6d 20 3d |..if (((|int)(m =|
|00001180| 20 28 73 74 72 75 63 74 | 20 6d 61 74 72 69 78 20 | (struct| matrix |
|00001190| 2a 29 6d 61 6c 6c 6f 63 | 28 73 69 7a 65 6f 66 28 |*)malloc|(sizeof(|
|000011a0| 73 74 72 75 63 74 20 6d | 61 74 72 69 78 29 20 2b |struct m|atrix) +|
|000011b0| 20 28 73 69 7a 65 6f 66 | 28 64 6f 75 62 6c 65 29 | (sizeof|(double)|
|000011c0| 20 2a 20 72 20 2a 20 63 | 29 29 29 29 20 3d 3d 20 | * r * c|)))) == |
|000011d0| 30 29 20 7b 5c 0d 0a 09 | 09 09 09 70 72 69 6e 74 |0) {\...|...print|
|000011e0| 66 28 22 41 6c 6c 6f 63 | 61 74 69 6f 6e 20 65 72 |f("Alloc|ation er|
|000011f0| 72 6f 72 3a 20 25 73 5c | 6e 22 2c 20 5f 5f 46 49 |ror: %s\|n", __FI|
|00001200| 4c 45 5f 5f 29 3b 5c 0d | 0a 09 09 09 09 65 78 69 |LE__);\.|.....exi|
|00001210| 74 28 31 29 3b 5c 0d 0a | 09 09 09 09 7d 5c 0d 0a |t(1);\..|....}\..|
|00001220| 09 09 09 6d 2d 3e 6d 5f | 72 6f 77 73 20 3d 20 72 |...m->m_|rows = r|
|00001230| 3b 5c 0d 0a 09 09 09 6d | 2d 3e 6d 5f 63 6f 6c 73 |;\.....m|->m_cols|
|00001240| 20 3d 20 63 3b 5c 0d 0a | 09 09 09 7d 0d 0a 2f 2f | = c;\..|...}..//|
|00001250| 67 6f 2e 73 79 73 69 6e | 20 64 64 20 2a 0d 0a 6d |go.sysin| dd *..m|
|00001260| 61 64 65 3d 54 52 55 45 | 0d 0a 69 66 20 5b 20 24 |ade=TRUE|..if [ $|
|00001270| 6d 61 64 65 20 3d 20 54 | 52 55 45 20 5d 3b 20 74 |made = T|RUE ]; t|
|00001280| 68 65 6e 0d 0a 09 2f 62 | 69 6e 2f 63 68 6d 6f 64 |hen.../b|in/chmod|
|00001290| 20 36 34 34 20 6d 61 74 | 2e 68 0d 0a 09 2f 62 69 | 644 mat|.h.../bi|
|000012a0| 6e 2f 65 63 68 6f 20 2d | 6e 20 27 09 27 3b 20 2f |n/echo -|n '.'; /|
|000012b0| 62 69 6e 2f 6c 73 20 2d | 6c 64 20 6d 61 74 2e 68 |bin/ls -|ld mat.h|
|000012c0| 0d 0a 66 69 0d 0a 2f 62 | 69 6e 2f 65 63 68 6f 20 |..fi../b|in/echo |
|000012d0| 27 45 78 74 72 61 63 74 | 69 6e 67 20 6d 5f 63 6f |'Extract|ing m_co|
|000012e0| 66 61 63 74 6f 72 2e 63 | 27 0d 0a 73 65 64 20 27 |factor.c|'..sed '|
|000012f0| 73 2f 5e 58 2f 2f 27 20 | 3c 3c 27 2f 2f 67 6f 2e |s/^X//' |<<'//go.|
|00001300| 73 79 73 69 6e 20 64 64 | 20 2a 27 20 3e 6d 5f 63 |sysin dd| *' >m_c|
|00001310| 6f 66 61 63 74 6f 72 2e | 63 0d 0a 73 74 61 74 69 |ofactor.|c..stati|
|00001320| 63 20 63 68 61 72 20 2a | 73 63 63 73 69 64 20 3d |c char *|sccsid =|
|00001330| 20 22 40 28 23 29 6d 5f | 63 6f 66 61 63 74 6f 72 | "@(#)m_|cofactor|
|00001340| 2e 63 09 34 2f 35 2f 38 | 32 20 28 55 20 6f 66 20 |.c.4/5/8|2 (U of |
|00001350| 4d 61 72 79 6c 61 6e 64 | 2c 20 46 4c 42 29 22 3b |Maryland|, FLB)";|
|00001360| 0d 0a 0d 0a 23 69 6e 63 | 6c 75 64 65 20 22 6d 61 |....#inc|lude "ma|
|00001370| 74 2e 68 22 0d 0a 0d 0a | 64 6f 75 62 6c 65 0d 0a |t.h"....|double..|
|00001380| 6d 5f 63 6f 66 61 63 74 | 6f 72 28 6d 61 74 2c 20 |m_cofact|or(mat, |
|00001390| 69 2c 20 6a 29 0d 0a 72 | 65 67 69 73 74 65 72 20 |i, j)..r|egister |
|000013a0| 73 74 72 75 63 74 20 6d | 61 74 72 69 78 20 2a 6d |struct m|atrix *m|
|000013b0| 61 74 3b 0d 0a 72 65 67 | 69 73 74 65 72 20 69 6e |at;..reg|ister in|
|000013c0| 74 20 69 2c 20 6a 3b 0d | 0a 7b 0d 0a 72 65 67 69 |t i, j;.|.{..regi|
|000013d0| 73 74 65 72 20 73 74 72 | 75 63 74 20 6d 61 74 72 |ster str|uct matr|
|000013e0| 69 78 20 2a 72 65 73 75 | 6c 74 3b 0d 0a 72 65 67 |ix *resu|lt;..reg|
|000013f0| 69 73 74 65 72 20 69 6e | 74 20 72 6f 77 2c 20 63 |ister in|t row, c|
|00001400| 6f 6c 2c 20 6f 5f 72 6f | 77 20 3d 20 30 2c 20 6f |ol, o_ro|w = 0, o|
|00001410| 5f 63 6f 6c 20 3d 20 30 | 3b 0d 0a 64 6f 75 62 6c |_col = 0|;..doubl|
|00001420| 65 20 64 65 74 3b 0d 0a | 0d 0a 6d 5f 63 72 65 61 |e det;..|..m_crea|
|00001430| 74 65 28 72 65 73 75 6c | 74 2c 20 6d 61 74 2d 3e |te(resul|t, mat->|
|00001440| 6d 5f 72 6f 77 73 20 2d | 20 31 2c 20 6d 61 74 2d |m_rows -| 1, mat-|
|00001450| 3e 6d 5f 63 6f 6c 73 20 | 2d 20 31 29 3b 0d 0a 0d |>m_cols |- 1);...|
|00001460| 0a 66 6f 72 20 28 72 6f | 77 20 3d 20 30 3b 20 72 |.for (ro|w = 0; r|
|00001470| 6f 77 20 3c 20 6d 61 74 | 2d 3e 6d 5f 72 6f 77 73 |ow < mat|->m_rows|
|00001480| 3b 20 72 6f 77 2b 2b 29 | 0d 0a 09 66 6f 72 20 28 |; row++)|...for (|
|00001490| 63 6f 6c 20 3d 20 30 3b | 20 63 6f 6c 20 3c 20 6d |col = 0;| col < m|
|000014a0| 61 74 2d 3e 6d 5f 63 6f | 6c 73 3b 20 63 6f 6c 2b |at->m_co|ls; col+|
|000014b0| 2b 29 0d 0a 09 09 69 66 | 20 28 72 6f 77 20 21 3d |+)....if| (row !=|
|000014c0| 20 69 20 26 26 20 63 6f | 6c 20 21 3d 20 6a 29 0d | i && co|l != j).|
|000014d0| 0a 09 09 09 6d 5f 76 28 | 72 65 73 75 6c 74 2c 20 |....m_v(|result, |
|000014e0| 6f 5f 72 6f 77 2b 2b 2c | 20 6f 5f 63 6f 6c 2b 2b |o_row++,| o_col++|
|000014f0| 29 20 3d 20 6d 5f 76 28 | 6d 61 74 2c 20 72 6f 77 |) = m_v(|mat, row|
|00001500| 2c 20 63 6f 6c 29 3b 0d | 0a 0d 0a 64 65 74 20 3d |, col);.|...det =|
|00001510| 20 6d 5f 64 65 74 65 72 | 6d 69 6e 61 6e 74 28 72 | m_deter|minant(r|
|00001520| 65 73 75 6c 74 29 3b 0d | 0a 66 72 65 65 28 72 65 |esult);.|.free(re|
|00001530| 73 75 6c 74 29 3b 0d 0a | 0d 0a 72 65 74 75 72 6e |sult);..|..return|
|00001540| 28 28 28 69 20 2b 20 6a | 29 20 26 20 30 31 29 3f |(((i + j|) & 01)?|
|00001550| 20 2d 64 65 74 3a 20 64 | 65 74 29 3b 0d 0a 7d 0d | -det: d|et);..}.|
|00001560| 0a 2f 2f 67 6f 2e 73 79 | 73 69 6e 20 64 64 20 2a |.//go.sy|sin dd *|
|00001570| 0d 0a 6d 61 64 65 3d 54 | 52 55 45 0d 0a 69 66 20 |..made=T|RUE..if |
|00001580| 5b 20 24 6d 61 64 65 20 | 3d 20 54 52 55 45 20 5d |[ $made |= TRUE ]|
|00001590| 3b 20 74 68 65 6e 0d 0a | 09 2f 62 69 6e 2f 63 68 |; then..|./bin/ch|
|000015a0| 6d 6f 64 20 36 34 34 20 | 6d 5f 63 6f 66 61 63 74 |mod 644 |m_cofact|
|000015b0| 6f 72 2e 63 0d 0a 09 2f | 62 69 6e 2f 65 63 68 6f |or.c.../|bin/echo|
|000015c0| 20 2d 6e 20 27 09 27 3b | 20 2f 62 69 6e 2f 6c 73 | -n '.';| /bin/ls|
|000015d0| 20 2d 6c 64 20 6d 5f 63 | 6f 66 61 63 74 6f 72 2e | -ld m_c|ofactor.|
|000015e0| 63 0d 0a 66 69 0d 0a 2f | 62 69 6e 2f 65 63 68 6f |c..fi../|bin/echo|
|000015f0| 20 27 45 78 74 72 61 63 | 74 69 6e 67 20 6d 5f 63 | 'Extrac|ting m_c|
|00001600| 6f 70 79 2e 63 27 0d 0a | 73 65 64 20 27 73 2f 5e |opy.c'..|sed 's/^|
|00001610| 58 2f 2f 27 20 3c 3c 27 | 2f 2f 67 6f 2e 73 79 73 |X//' <<'|//go.sys|
|00001620| 69 6e 20 64 64 20 2a 27 | 20 3e 6d 5f 63 6f 70 79 |in dd *'| >m_copy|
|00001630| 2e 63 0d 0a 73 74 61 74 | 69 63 20 63 68 61 72 20 |.c..stat|ic char |
|00001640| 2a 73 63 63 73 69 64 20 | 3d 20 22 40 28 23 29 6d |*sccsid |= "@(#)m|
|00001650| 5f 63 6f 70 79 2e 63 09 | 34 2f 35 2f 38 32 20 28 |_copy.c.|4/5/82 (|
|00001660| 55 20 6f 66 20 4d 61 72 | 79 6c 61 6e 64 2c 20 46 |U of Mar|yland, F|
|00001670| 4c 42 29 22 3b 0d 0a 0d | 0a 23 69 6e 63 6c 75 64 |LB)";...|.#includ|
|00001680| 65 20 22 6d 61 74 2e 68 | 22 0d 0a 0d 0a 73 74 72 |e "mat.h|"....str|
|00001690| 75 63 74 20 6d 61 74 72 | 69 78 20 2a 0d 0a 6d 5f |uct matr|ix *..m_|
|000016a0| 63 6f 70 79 28 6d 61 74 | 29 0d 0a 72 65 67 69 73 |copy(mat|)..regis|
|000016b0| 74 65 72 20 73 74 72 75 | 63 74 20 6d 61 74 72 69 |ter stru|ct matri|
|000016c0| 78 20 2a 6d 61 74 3b 0d | 0a 7b 0d 0a 72 65 67 69 |x *mat;.|.{..regi|
|000016d0| 73 74 65 72 20 73 74 72 | 75 63 74 20 6d 61 74 72 |ster str|uct matr|
|000016e0| 69 78 20 2a 72 65 73 75 | 6c 74 3b 0d 0a 72 65 67 |ix *resu|lt;..reg|
|000016f0| 69 73 74 65 72 20 69 6e | 74 20 72 6f 77 2c 20 63 |ister in|t row, c|
|00001700| 6f 6c 3b 0d 0a 0d 0a 6d | 5f 63 72 65 61 74 65 28 |ol;....m|_create(|
|00001710| 72 65 73 75 6c 74 2c 20 | 6d 61 74 2d 3e 6d 5f 72 |result, |mat->m_r|
|00001720| 6f 77 73 2c 20 6d 61 74 | 2d 3e 6d 5f 63 6f 6c 73 |ows, mat|->m_cols|
|00001730| 29 3b 0d 0a 0d 0a 66 6f | 72 20 28 72 6f 77 20 3d |);....fo|r (row =|
|00001740| 20 30 3b 20 72 6f 77 20 | 3c 20 6d 61 74 2d 3e 6d | 0; row |< mat->m|
|00001750| 5f 72 6f 77 73 3b 20 72 | 6f 77 2b 2b 29 0d 0a 09 |_rows; r|ow++)...|
|00001760| 66 6f 72 20 28 63 6f 6c | 20 3d 20 30 3b 20 63 6f |for (col| = 0; co|
|00001770| 6c 20 3c 20 6d 61 74 2d | 3e 6d 5f 63 6f 6c 73 3b |l < mat-|>m_cols;|
|00001780| 20 63 6f 6c 2b 2b 29 0d | 0a 09 09 6d 5f 76 28 72 | col++).|...m_v(r|
|00001790| 65 73 75 6c 74 2c 20 72 | 6f 77 2c 20 63 6f 6c 29 |esult, r|ow, col)|
|000017a0| 20 3d 20 6d 5f 76 28 6d | 61 74 2c 20 72 6f 77 2c | = m_v(m|at, row,|
|000017b0| 20 63 6f 6c 29 3b 0d 0a | 0d 0a 72 65 74 75 72 6e | col);..|..return|
|000017c0| 28 72 65 73 75 6c 74 29 | 3b 0d 0a 7d 0d 0a 2f 2f |(result)|;..}..//|
|000017d0| 67 6f 2e 73 79 73 69 6e | 20 64 64 20 2a 0d 0a 6d |go.sysin| dd *..m|
|000017e0| 61 64 65 3d 54 52 55 45 | 0d 0a 69 66 20 5b 20 24 |ade=TRUE|..if [ $|
|000017f0| 6d 61 64 65 20 3d 20 54 | 52 55 45 20 5d 3b 20 74 |made = T|RUE ]; t|
|00001800| 68 65 6e 0d 0a 09 2f 62 | 69 6e 2f 63 68 6d 6f 64 |hen.../b|in/chmod|
|00001810| 20 36 34 34 20 6d 5f 63 | 6f 70 79 2e 63 0d 0a 09 | 644 m_c|opy.c...|
|00001820| 2f 62 69 6e 2f 65 63 68 | 6f 20 2d 6e 20 27 09 27 |/bin/ech|o -n '.'|
|00001830| 3b 20 2f 62 69 6e 2f 6c | 73 20 2d 6c 64 20 6d 5f |; /bin/l|s -ld m_|
|00001840| 63 6f 70 79 2e 63 0d 0a | 66 69 0d 0a 2f 62 69 6e |copy.c..|fi../bin|
|00001850| 2f 65 63 68 6f 20 27 45 | 78 74 72 61 63 74 69 6e |/echo 'E|xtractin|
|00001860| 67 20 6d 5f 64 65 74 65 | 72 6d 2e 63 27 0d 0a 73 |g m_dete|rm.c'..s|
|00001870| 65 64 20 27 73 2f 5e 58 | 2f 2f 27 20 3c 3c 27 2f |ed 's/^X|//' <<'/|
|00001880| 2f 67 6f 2e 73 79 73 69 | 6e 20 64 64 20 2a 27 20 |/go.sysi|n dd *' |
|00001890| 3e 6d 5f 64 65 74 65 72 | 6d 2e 63 0d 0a 73 74 61 |>m_deter|m.c..sta|
|000018a0| 74 69 63 20 63 68 61 72 | 20 2a 73 63 63 73 69 64 |tic char| *sccsid|
|000018b0| 20 3d 20 22 40 28 23 29 | 6d 5f 64 65 74 65 72 6d | = "@(#)|m_determ|
|000018c0| 69 6e 61 6e 74 2e 63 09 | 34 2f 35 2f 38 32 20 28 |inant.c.|4/5/82 (|
|000018d0| 55 20 6f 66 20 4d 61 72 | 79 6c 61 6e 64 2c 20 46 |U of Mar|yland, F|
|000018e0| 4c 42 29 22 3b 0d 0a 0d | 0a 23 69 6e 63 6c 75 64 |LB)";...|.#includ|
|000018f0| 65 20 22 6d 61 74 2e 68 | 22 0d 0a 0d 0a 64 6f 75 |e "mat.h|"....dou|
|00001900| 62 6c 65 0d 0a 6d 5f 64 | 65 74 65 72 6d 69 6e 61 |ble..m_d|etermina|
|00001910| 6e 74 28 6d 61 74 29 0d | 0a 72 65 67 69 73 74 65 |nt(mat).|.registe|
|00001920| 72 20 73 74 72 75 63 74 | 20 6d 61 74 72 69 78 20 |r struct| matrix |
|00001930| 2a 6d 61 74 3b 0d 0a 7b | 0d 0a 72 65 67 69 73 74 |*mat;..{|..regist|
|00001940| 65 72 20 69 6e 74 20 63 | 6f 6c 3b 0d 0a 64 6f 75 |er int c|ol;..dou|
|00001950| 62 6c 65 20 64 65 74 20 | 3d 20 30 2e 30 3b 0d 0a |ble det |= 0.0;..|
|00001960| 0d 0a 69 66 20 28 6d 61 | 74 2d 3e 6d 5f 72 6f 77 |..if (ma|t->m_row|
|00001970| 73 20 3d 3d 20 31 29 0d | 0a 09 72 65 74 75 72 6e |s == 1).|..return|
|00001980| 28 6d 5f 76 28 6d 61 74 | 2c 20 30 2c 20 30 29 29 |(m_v(mat|, 0, 0))|
|00001990| 3b 0d 0a 0d 0a 66 6f 72 | 20 28 63 6f 6c 20 3d 20 |;....for| (col = |
|000019a0| 30 3b 20 63 6f 6c 20 3c | 20 6d 61 74 2d 3e 6d 5f |0; col <| mat->m_|
|000019b0| 63 6f 6c 73 3b 20 63 6f | 6c 2b 2b 29 0d 0a 09 64 |cols; co|l++)...d|
|000019c0| 65 74 20 2b 3d 20 09 6d | 5f 76 28 6d 61 74 2c 20 |et += .m|_v(mat, |
|000019d0| 30 2c 20 63 6f 6c 29 20 | 2a 20 6d 5f 63 6f 66 61 |0, col) |* m_cofa|
|000019e0| 63 74 6f 72 28 6d 61 74 | 2c 20 30 2c 20 63 6f 6c |ctor(mat|, 0, col|
|000019f0| 29 3b 0d 0a 0d 0a 72 65 | 74 75 72 6e 28 64 65 74 |);....re|turn(det|
|00001a00| 29 3b 0d 0a 7d 0d 0a 2f | 2f 67 6f 2e 73 79 73 69 |);..}../|/go.sysi|
|00001a10| 6e 20 64 64 20 2a 0d 0a | 6d 61 64 65 3d 54 52 55 |n dd *..|made=TRU|
|00001a20| 45 0d 0a 69 66 20 5b 20 | 24 6d 61 64 65 20 3d 20 |E..if [ |$made = |
|00001a30| 54 52 55 45 20 5d 3b 20 | 74 68 65 6e 0d 0a 09 2f |TRUE ]; |then.../|
|00001a40| 62 69 6e 2f 63 68 6d 6f | 64 20 36 34 34 20 6d 5f |bin/chmo|d 644 m_|
|00001a50| 64 65 74 65 72 6d 2e 63 | 0d 0a 09 2f 62 69 6e 2f |determ.c|.../bin/|
|00001a60| 65 63 68 6f 20 2d 6e 20 | 27 09 27 3b 20 2f 62 69 |echo -n |'.'; /bi|
|00001a70| 6e 2f 6c 73 20 2d 6c 64 | 20 6d 5f 64 65 74 65 72 |n/ls -ld| m_deter|
|00001a80| 6d 2e 63 0d 0a 66 69 0d | 0a 2f 62 69 6e 2f 65 63 |m.c..fi.|./bin/ec|
|00001a90| 68 6f 20 27 45 78 74 72 | 61 63 74 69 6e 67 20 6d |ho 'Extr|acting m|
|00001aa0| 5f 64 75 6d 70 2e 63 27 | 0d 0a 73 65 64 20 27 73 |_dump.c'|..sed 's|
|00001ab0| 2f 5e 58 2f 2f 27 20 3c | 3c 27 2f 2f 67 6f 2e 73 |/^X//' <|<'//go.s|
|00001ac0| 79 73 69 6e 20 64 64 20 | 2a 27 20 3e 6d 5f 64 75 |ysin dd |*' >m_du|
|00001ad0| 6d 70 2e 63 0d 0a 73 74 | 61 74 69 63 20 63 68 61 |mp.c..st|atic cha|
|00001ae0| 72 20 2a 73 63 63 73 69 | 64 20 3d 20 22 40 28 23 |r *sccsi|d = "@(#|
|00001af0| 29 6d 5f 64 75 6d 70 2e | 63 09 34 2f 36 2f 38 32 |)m_dump.|c.4/6/82|
|00001b00| 20 28 55 20 6f 66 20 4d | 61 72 79 6c 61 6e 64 2c | (U of M|aryland,|
|00001b10| 20 46 4c 42 29 22 3b 0d | 0a 0d 0a 23 69 6e 63 6c | FLB)";.|...#incl|
|00001b20| 75 64 65 20 22 6d 61 74 | 2e 68 22 0d 0a 0d 0a 73 |ude "mat|.h"....s|
|00001b30| 74 72 75 63 74 20 6d 61 | 74 72 69 78 20 2a 0d 0a |truct ma|trix *..|
|00001b40| 6d 5f 64 75 6d 70 28 6d | 61 74 29 0d 0a 72 65 67 |m_dump(m|at)..reg|
|00001b50| 69 73 74 65 72 20 73 74 | 72 75 63 74 20 6d 61 74 |ister st|ruct mat|
|00001b60| 72 69 78 20 2a 6d 61 74 | 3b 0d 0a 7b 0d 0a 72 65 |rix *mat|;..{..re|
|00001b70| 67 69 73 74 65 72 20 69 | 6e 74 20 72 6f 77 2c 20 |gister i|nt row, |
|00001b80| 63 6f 6c 3b 0d 0a 0d 0a | 70 72 69 6e 74 66 28 22 |col;....|printf("|
|00001b90| 25 64 20 58 20 25 64 5c | 6e 22 2c 20 6d 61 74 2d |%d X %d\|n", mat-|
|00001ba0| 3e 6d 5f 72 6f 77 73 2c | 20 6d 61 74 2d 3e 6d 5f |>m_rows,| mat->m_|
|00001bb0| 63 6f 6c 73 29 3b 0d 0a | 0d 0a 66 6f 72 20 28 72 |cols);..|..for (r|
|00001bc0| 6f 77 20 3d 20 30 3b 20 | 72 6f 77 20 3c 20 6d 61 |ow = 0; |row < ma|
|00001bd0| 74 2d 3e 6d 5f 72 6f 77 | 73 3b 20 72 6f 77 2b 2b |t->m_row|s; row++|
|00001be0| 29 20 7b 0d 0a 09 66 6f | 72 20 28 63 6f 6c 20 3d |) {...fo|r (col =|
|00001bf0| 20 30 3b 20 63 6f 6c 20 | 3c 20 6d 61 74 2d 3e 6d | 0; col |< mat->m|
|00001c00| 5f 63 6f 6c 73 3b 20 63 | 6f 6c 2b 2b 29 0d 0a 09 |_cols; c|ol++)...|
|00001c10| 09 70 72 69 6e 74 66 28 | 22 25 66 2c 20 22 2c 20 |.printf(|"%f, ", |
|00001c20| 6d 5f 76 28 6d 61 74 2c | 20 72 6f 77 2c 20 63 6f |m_v(mat,| row, co|
|00001c30| 6c 29 29 3b 0d 0a 09 70 | 75 74 63 68 61 72 28 27 |l));...p|utchar('|
|00001c40| 5c 6e 27 29 3b 0d 0a 09 | 7d 0d 0a 7d 0d 0a 2f 2f |\n');...|}..}..//|
|00001c50| 67 6f 2e 73 79 73 69 6e | 20 64 64 20 2a 0d 0a 6d |go.sysin| dd *..m|
|00001c60| 61 64 65 3d 54 52 55 45 | 0d 0a 69 66 20 5b 20 24 |ade=TRUE|..if [ $|
|00001c70| 6d 61 64 65 20 3d 20 54 | 52 55 45 20 5d 3b 20 74 |made = T|RUE ]; t|
|00001c80| 68 65 6e 0d 0a 09 2f 62 | 69 6e 2f 63 68 6d 6f 64 |hen.../b|in/chmod|
|00001c90| 20 36 34 34 20 6d 5f 64 | 75 6d 70 2e 63 0d 0a 09 | 644 m_d|ump.c...|
|00001ca0| 2f 62 69 6e 2f 65 63 68 | 6f 20 2d 6e 20 27 09 27 |/bin/ech|o -n '.'|
|00001cb0| 3b 20 2f 62 69 6e 2f 6c | 73 20 2d 6c 64 20 6d 5f |; /bin/l|s -ld m_|
|00001cc0| 64 75 6d 70 2e 63 0d 0a | 66 69 0d 0a 2f 62 69 6e |dump.c..|fi../bin|
|00001cd0| 2f 65 63 68 6f 20 27 45 | 78 74 72 61 63 74 69 6e |/echo 'E|xtractin|
|00001ce0| 67 20 6d 5f 69 6e 76 65 | 72 74 2e 63 27 0d 0a 73 |g m_inve|rt.c'..s|
|00001cf0| 65 64 20 27 73 2f 5e 58 | 2f 2f 27 20 3c 3c 27 2f |ed 's/^X|//' <<'/|
|00001d00| 2f 67 6f 2e 73 79 73 69 | 6e 20 64 64 20 2a 27 20 |/go.sysi|n dd *' |
|00001d10| 3e 6d 5f 69 6e 76 65 72 | 74 2e 63 0d 0a 73 74 61 |>m_inver|t.c..sta|
|00001d20| 74 69 63 20 63 68 61 72 | 20 2a 73 63 63 73 69 64 |tic char| *sccsid|
|00001d30| 20 3d 20 22 40 28 23 29 | 6d 5f 69 6e 76 65 72 74 | = "@(#)|m_invert|
|00001d40| 2e 63 09 34 2f 37 2f 38 | 32 20 28 55 20 6f 66 20 |.c.4/7/8|2 (U of |
|00001d50| 4d 61 72 79 6c 61 6e 64 | 2c 20 46 4c 42 29 22 3b |Maryland|, FLB)";|
|00001d60| 0d 0a 0d 0a 23 69 6e 63 | 6c 75 64 65 20 22 6d 61 |....#inc|lude "ma|
|00001d70| 74 2e 68 22 0d 0a 0d 0a | 73 74 72 75 63 74 20 6d |t.h"....|struct m|
|00001d80| 61 74 72 69 78 20 2a 0d | 0a 6d 5f 69 6e 76 65 72 |atrix *.|.m_inver|
|00001d90| 74 28 6d 61 74 29 0d 0a | 72 65 67 69 73 74 65 72 |t(mat)..|register|
|00001da0| 20 73 74 72 75 63 74 20 | 6d 61 74 72 69 78 20 2a | struct |matrix *|
|00001db0| 6d 61 74 3b 0d 0a 7b 0d | 0a 72 65 67 69 73 74 65 |mat;..{.|.registe|
|00001dc0| 72 20 73 74 72 75 63 74 | 20 6d 61 74 72 69 78 20 |r struct| matrix |
|00001dd0| 2a 72 65 73 75 6c 74 3b | 0d 0a 72 65 67 69 73 74 |*result;|..regist|
|00001de0| 65 72 20 69 6e 74 20 72 | 6f 77 2c 20 63 6f 6c 3b |er int r|ow, col;|
|00001df0| 0d 0a 64 6f 75 62 6c 65 | 20 64 65 74 3b 0d 0a 0d |..double| det;...|
|00001e00| 0a 69 66 28 28 64 65 74 | 20 3d 20 6d 5f 64 65 74 |.if((det| = m_det|
|00001e10| 65 72 6d 69 6e 61 6e 74 | 28 6d 61 74 29 29 20 3d |erminant|(mat)) =|
|00001e20| 3d 20 30 2e 30 29 20 7b | 0d 0a 09 70 72 69 6e 74 |= 0.0) {|...print|
|00001e30| 66 28 22 6d 5f 69 6e 76 | 65 72 74 3a 20 73 69 6e |f("m_inv|ert: sin|
|00001e40| 67 75 6c 61 72 20 6d 61 | 74 72 69 78 5c 6e 22 29 |gular ma|trix\n")|
|00001e50| 3b 0d 0a 09 72 65 74 75 | 72 6e 28 4d 5f 4e 55 4c |;...retu|rn(M_NUL|
|00001e60| 4c 29 3b 0d 0a 09 7d 0d | 0a 0d 0a 6d 5f 63 72 65 |L);...}.|...m_cre|
|00001e70| 61 74 65 28 72 65 73 75 | 6c 74 2c 20 6d 61 74 2d |ate(resu|lt, mat-|
|00001e80| 3e 6d 5f 63 6f 6c 73 2c | 20 6d 61 74 2d 3e 6d 5f |>m_cols,| mat->m_|
|00001e90| 72 6f 77 73 29 3b 0d 0a | 0d 0a 66 6f 72 20 28 72 |rows);..|..for (r|
|00001ea0| 6f 77 20 3d 20 30 3b 20 | 72 6f 77 20 3c 20 6d 61 |ow = 0; |row < ma|
|00001eb0| 74 2d 3e 6d 5f 72 6f 77 | 73 3b 20 72 6f 77 2b 2b |t->m_row|s; row++|
|00001ec0| 29 0d 0a 09 66 6f 72 20 | 28 63 6f 6c 20 3d 20 30 |)...for |(col = 0|
|00001ed0| 3b 20 63 6f 6c 20 3c 20 | 6d 61 74 2d 3e 6d 5f 63 |; col < |mat->m_c|
|00001ee0| 6f 6c 73 3b 20 63 6f 6c | 2b 2b 29 0d 0a 09 09 6d |ols; col|++)....m|
|00001ef0| 5f 76 28 72 65 73 75 6c | 74 2c 20 63 6f 6c 2c 20 |_v(resul|t, col, |
|00001f00| 72 6f 77 29 20 3d 20 6d | 5f 63 6f 66 61 63 74 6f |row) = m|_cofacto|
|00001f10| 72 28 6d 61 74 2c 20 72 | 6f 77 2c 20 63 6f 6c 29 |r(mat, r|ow, col)|
|00001f20| 20 2f 20 64 65 74 3b 0d | 0a 0d 0a 72 65 74 75 72 | / det;.|...retur|
|00001f30| 6e 28 72 65 73 75 6c 74 | 29 3b 0d 0a 7d 0d 0a 2f |n(result|);..}../|
|00001f40| 2f 67 6f 2e 73 79 73 69 | 6e 20 64 64 20 2a 0d 0a |/go.sysi|n dd *..|
|00001f50| 6d 61 64 65 3d 54 52 55 | 45 0d 0a 69 66 20 5b 20 |made=TRU|E..if [ |
|00001f60| 24 6d 61 64 65 20 3d 20 | 54 52 55 45 20 5d 3b 20 |$made = |TRUE ]; |
|00001f70| 74 68 65 6e 0d 0a 09 2f | 62 69 6e 2f 63 68 6d 6f |then.../|bin/chmo|
|00001f80| 64 20 36 34 34 20 6d 5f | 69 6e 76 65 72 74 2e 63 |d 644 m_|invert.c|
|00001f90| 0d 0a 09 2f 62 69 6e 2f | 65 63 68 6f 20 2d 6e 20 |.../bin/|echo -n |
|00001fa0| 27 09 27 3b 20 2f 62 69 | 6e 2f 6c 73 20 2d 6c 64 |'.'; /bi|n/ls -ld|
|00001fb0| 20 6d 5f 69 6e 76 65 72 | 74 2e 63 0d 0a 66 69 0d | m_inver|t.c..fi.|
|00001fc0| 0a 2f 62 69 6e 2f 65 63 | 68 6f 20 27 45 78 74 72 |./bin/ec|ho 'Extr|
|00001fd0| 61 63 74 69 6e 67 20 6d | 5f 6d 75 6c 74 69 70 6c |acting m|_multipl|
|00001fe0| 79 2e 63 27 0d 0a 73 65 | 64 20 27 73 2f 5e 58 2f |y.c'..se|d 's/^X/|
|00001ff0| 2f 27 20 3c 3c 27 2f 2f | 67 6f 2e 73 79 73 69 6e |/' <<'//|go.sysin|
|00002000| 20 64 64 20 2a 27 20 3e | 6d 5f 6d 75 6c 74 69 70 | dd *' >|m_multip|
|00002010| 6c 79 2e 63 0d 0a 73 74 | 61 74 69 63 20 63 68 61 |ly.c..st|atic cha|
|00002020| 72 20 2a 73 63 63 73 69 | 64 20 3d 20 22 40 28 23 |r *sccsi|d = "@(#|
|00002030| 29 6d 5f 6d 75 6c 74 69 | 70 6c 79 2e 63 09 34 2f |)m_multi|ply.c.4/|
|00002040| 36 2f 38 32 20 28 55 20 | 6f 66 20 4d 61 72 79 6c |6/82 (U |of Maryl|
|00002050| 61 6e 64 2c 20 46 4c 42 | 29 22 3b 0d 0a 0d 0a 23 |and, FLB|)";....#|
|00002060| 69 6e 63 6c 75 64 65 20 | 22 6d 61 74 2e 68 22 0d |include |"mat.h".|
|00002070| 0a 0d 0a 73 74 72 75 63 | 74 20 6d 61 74 72 69 78 |...struc|t matrix|
|00002080| 20 2a 0d 0a 6d 5f 6d 75 | 6c 74 69 70 6c 79 28 6d | *..m_mu|ltiply(m|
|00002090| 61 74 31 2c 20 6d 61 74 | 32 29 0d 0a 72 65 67 69 |at1, mat|2)..regi|
|000020a0| 73 74 65 72 20 73 74 72 | 75 63 74 20 6d 61 74 72 |ster str|uct matr|
|000020b0| 69 78 20 2a 6d 61 74 31 | 2c 20 2a 6d 61 74 32 3b |ix *mat1|, *mat2;|
|000020c0| 0d 0a 7b 0d 0a 72 65 67 | 69 73 74 65 72 20 73 74 |..{..reg|ister st|
|000020d0| 72 75 63 74 20 6d 61 74 | 72 69 78 20 2a 72 65 73 |ruct mat|rix *res|
|000020e0| 75 6c 74 3b 0d 0a 72 65 | 67 69 73 74 65 72 20 69 |ult;..re|gister i|
|000020f0| 6e 74 20 72 6f 77 2c 20 | 63 6f 6c 2c 20 69 78 3b |nt row, |col, ix;|
|00002100| 0d 0a 64 6f 75 62 6c 65 | 20 73 75 6d 3b 0d 0a 0d |..double| sum;...|
|00002110| 0a 69 66 20 28 6d 61 74 | 31 2d 3e 6d 5f 63 6f 6c |.if (mat|1->m_col|
|00002120| 73 20 21 3d 20 6d 61 74 | 32 2d 3e 6d 5f 72 6f 77 |s != mat|2->m_row|
|00002130| 73 29 20 7b 0d 0a 09 70 | 72 69 6e 74 66 28 22 6d |s) {...p|rintf("m|
|00002140| 5f 6d 75 6c 74 69 70 6c | 79 3a 20 6d 61 74 72 69 |_multipl|y: matri|
|00002150| 63 65 73 20 6e 6f 74 20 | 63 6f 6d 70 61 74 69 62 |ces not |compatib|
|00002160| 6c 65 2e 5c 6e 22 29 3b | 0d 0a 09 72 65 74 75 72 |le.\n");|...retur|
|00002170| 6e 28 4d 5f 4e 55 4c 4c | 29 3b 0d 0a 09 7d 0d 0a |n(M_NULL|);...}..|
|00002180| 0d 0a 6d 5f 63 72 65 61 | 74 65 28 72 65 73 75 6c |..m_crea|te(resul|
|00002190| 74 2c 20 6d 61 74 31 2d | 3e 6d 5f 72 6f 77 73 2c |t, mat1-|>m_rows,|
|000021a0| 20 6d 61 74 32 2d 3e 6d | 5f 63 6f 6c 73 29 3b 0d | mat2->m|_cols);.|
|000021b0| 0a 0d 0a 66 6f 72 20 28 | 72 6f 77 20 3d 20 30 3b |...for (|row = 0;|
|000021c0| 20 72 6f 77 20 3c 20 6d | 61 74 31 2d 3e 6d 5f 72 | row < m|at1->m_r|
|000021d0| 6f 77 73 3b 20 72 6f 77 | 2b 2b 29 0d 0a 09 66 6f |ows; row|++)...fo|
|000021e0| 72 20 28 63 6f 6c 20 3d | 20 30 3b 20 63 6f 6c 20 |r (col =| 0; col |
|000021f0| 3c 20 6d 61 74 32 2d 3e | 6d 5f 63 6f 6c 73 3b 20 |< mat2->|m_cols; |
|00002200| 63 6f 6c 2b 2b 29 20 7b | 0d 0a 09 09 73 75 6d 20 |col++) {|....sum |
|00002210| 3d 20 30 2e 30 3b 0d 0a | 09 09 66 6f 72 20 28 69 |= 0.0;..|..for (i|
|00002220| 78 20 3d 20 30 3b 20 69 | 78 20 3c 20 6d 61 74 31 |x = 0; i|x < mat1|
|00002230| 2d 3e 6d 5f 63 6f 6c 73 | 3b 20 69 78 2b 2b 29 0d |->m_cols|; ix++).|
|00002240| 0a 09 09 09 73 75 6d 20 | 2b 3d 20 6d 5f 76 28 6d |....sum |+= m_v(m|
|00002250| 61 74 31 2c 20 72 6f 77 | 2c 20 69 78 29 20 2a 20 |at1, row|, ix) * |
|00002260| 6d 5f 76 28 6d 61 74 32 | 2c 20 69 78 2c 20 63 6f |m_v(mat2|, ix, co|
|00002270| 6c 29 3b 0d 0a 09 09 6d | 5f 76 28 72 65 73 75 6c |l);....m|_v(resul|
|00002280| 74 2c 20 72 6f 77 2c 20 | 63 6f 6c 29 20 3d 20 73 |t, row, |col) = s|
|00002290| 75 6d 3b 0d 0a 09 09 7d | 0d 0a 0d 0a 72 65 74 75 |um;....}|....retu|
|000022a0| 72 6e 28 72 65 73 75 6c | 74 29 3b 0d 0a 7d 0d 0a |rn(resul|t);..}..|
|000022b0| 2f 2f 67 6f 2e 73 79 73 | 69 6e 20 64 64 20 2a 0d |//go.sys|in dd *.|
|000022c0| 0a 6d 61 64 65 3d 54 52 | 55 45 0d 0a 69 66 20 5b |.made=TR|UE..if [|
|000022d0| 20 24 6d 61 64 65 20 3d | 20 54 52 55 45 20 5d 3b | $made =| TRUE ];|
|000022e0| 20 74 68 65 6e 0d 0a 09 | 2f 62 69 6e 2f 63 68 6d | then...|/bin/chm|
|000022f0| 6f 64 20 36 34 34 20 6d | 5f 6d 75 6c 74 69 70 6c |od 644 m|_multipl|
|00002300| 79 2e 63 0d 0a 09 2f 62 | 69 6e 2f 65 63 68 6f 20 |y.c.../b|in/echo |
|00002310| 2d 6e 20 27 09 27 3b 20 | 2f 62 69 6e 2f 6c 73 20 |-n '.'; |/bin/ls |
|00002320| 2d 6c 64 20 6d 5f 6d 75 | 6c 74 69 70 6c 79 2e 63 |-ld m_mu|ltiply.c|
|00002330| 0d 0a 66 69 0d 0a 2f 62 | 69 6e 2f 65 63 68 6f 20 |..fi../b|in/echo |
|00002340| 27 45 78 74 72 61 63 74 | 69 6e 67 20 6d 5f 72 65 |'Extract|ing m_re|
|00002350| 61 64 2e 63 27 0d 0a 73 | 65 64 20 27 73 2f 5e 58 |ad.c'..s|ed 's/^X|
|00002360| 2f 2f 27 20 3c 3c 27 2f | 2f 67 6f 2e 73 79 73 69 |//' <<'/|/go.sysi|
|00002370| 6e 20 64 64 20 2a 27 20 | 3e 6d 5f 72 65 61 64 2e |n dd *' |>m_read.|
|00002380| 63 0d 0a 73 74 61 74 69 | 63 20 63 68 61 72 20 2a |c..stati|c char *|
|00002390| 73 63 63 73 69 64 20 3d | 20 22 40 28 23 29 6d 5f |sccsid =| "@(#)m_|
|000023a0| 72 65 61 64 2e 63 09 34 | 2f 36 2f 38 32 20 28 55 |read.c.4|/6/82 (U|
|000023b0| 20 6f 66 20 4d 61 72 79 | 6c 61 6e 64 2c 20 46 4c | of Mary|land, FL|
|000023c0| 42 29 22 3b 0d 0a 0d 0a | 23 69 6e 63 6c 75 64 65 |B)";....|#include|
|000023d0| 20 3c 73 74 64 69 6f 2e | 68 3e 0d 0a 23 69 6e 63 | <stdio.|h>..#inc|
|000023e0| 6c 75 64 65 20 22 6d 61 | 74 2e 68 22 0d 0a 0d 0a |lude "ma|t.h"....|
|000023f0| 73 74 72 75 63 74 20 6d | 61 74 72 69 78 20 2a 0d |struct m|atrix *.|
|00002400| 0a 6d 5f 72 65 61 64 28 | 29 0d 0a 7b 0d 0a 72 65 |.m_read(|)..{..re|
|00002410| 67 69 73 74 65 72 20 73 | 74 72 75 63 74 20 6d 61 |gister s|truct ma|
|00002420| 74 72 69 78 20 2a 72 65 | 73 75 6c 74 3b 0d 0a 72 |trix *re|sult;..r|
|00002430| 65 67 69 73 74 65 72 20 | 69 6e 74 20 72 6f 77 2c |egister |int row,|
|00002440| 20 63 6f 6c 3b 0d 0a 69 | 6e 74 20 72 6f 77 73 2c | col;..i|nt rows,|
|00002450| 20 63 6f 6c 73 3b 0d 0a | 0d 0a 73 63 61 6e 66 28 | cols;..|..scanf(|
|00002460| 22 25 64 25 64 22 2c 20 | 26 72 6f 77 73 2c 20 26 |"%d%d", |&rows, &|
|00002470| 63 6f 6c 73 29 3b 0d 0a | 6d 5f 63 72 65 61 74 65 |cols);..|m_create|
|00002480| 28 72 65 73 75 6c 74 2c | 20 72 6f 77 73 2c 20 63 |(result,| rows, c|
|00002490| 6f 6c 73 29 3b 0d 0a 0d | 0a 66 6f 72 20 28 72 6f |ols);...|.for (ro|
|000024a0| 77 20 3d 20 30 3b 20 72 | 6f 77 20 3c 20 72 6f 77 |w = 0; r|ow < row|
|000024b0| 73 3b 20 72 6f 77 2b 2b | 29 0d 0a 09 66 6f 72 20 |s; row++|)...for |
|000024c0| 28 63 6f 6c 20 3d 20 30 | 3b 20 63 6f 6c 20 3c 20 |(col = 0|; col < |
|000024d0| 63 6f 6c 73 3b 20 63 6f | 6c 2b 2b 29 0d 0a 09 09 |cols; co|l++)....|
|000024e0| 73 63 61 6e 66 28 22 25 | 6c 66 22 2c 20 26 6d 5f |scanf("%|lf", &m_|
|000024f0| 76 28 72 65 73 75 6c 74 | 2c 20 72 6f 77 2c 20 63 |v(result|, row, c|
|00002500| 6f 6c 29 29 3b 0d 0a 0d | 0a 72 65 74 75 72 6e 28 |ol));...|.return(|
|00002510| 72 65 73 75 6c 74 29 3b | 0d 0a 7d 0d 0a 2f 2f 67 |result);|..}..//g|
|00002520| 6f 2e 73 79 73 69 6e 20 | 64 64 20 2a 0d 0a 6d 61 |o.sysin |dd *..ma|
|00002530| 64 65 3d 54 52 55 45 0d | 0a 69 66 20 5b 20 24 6d |de=TRUE.|.if [ $m|
|00002540| 61 64 65 20 3d 20 54 52 | 55 45 20 5d 3b 20 74 68 |ade = TR|UE ]; th|
|00002550| 65 6e 0d 0a 09 2f 62 69 | 6e 2f 63 68 6d 6f 64 20 |en.../bi|n/chmod |
|00002560| 36 34 34 20 6d 5f 72 65 | 61 64 2e 63 0d 0a 09 2f |644 m_re|ad.c.../|
|00002570| 62 69 6e 2f 65 63 68 6f | 20 2d 6e 20 27 09 27 3b |bin/echo| -n '.';|
|00002580| 20 2f 62 69 6e 2f 6c 73 | 20 2d 6c 64 20 6d 5f 72 | /bin/ls| -ld m_r|
|00002590| 65 61 64 2e 63 0d 0a 66 | 69 0d 0a 2f 62 69 6e 2f |ead.c..f|i../bin/|
|000025a0| 65 63 68 6f 20 27 45 78 | 74 72 61 63 74 69 6e 67 |echo 'Ex|tracting|
|000025b0| 20 6d 5f 73 6f 6c 76 65 | 2e 63 27 0d 0a 73 65 64 | m_solve|.c'..sed|
|000025c0| 20 27 73 2f 5e 58 2f 2f | 27 20 3c 3c 27 2f 2f 67 | 's/^X//|' <<'//g|
|000025d0| 6f 2e 73 79 73 69 6e 20 | 64 64 20 2a 27 20 3e 6d |o.sysin |dd *' >m|
|000025e0| 5f 73 6f 6c 76 65 2e 63 | 0d 0a 73 74 61 74 69 63 |_solve.c|..static|
|000025f0| 20 63 68 61 72 20 2a 73 | 63 63 73 69 64 20 3d 20 | char *s|ccsid = |
|00002600| 22 40 28 23 29 6d 5f 73 | 6f 6c 76 65 2e 63 09 34 |"@(#)m_s|olve.c.4|
|00002610| 2f 31 36 2f 38 32 20 28 | 55 20 6f 66 20 4d 61 72 |/16/82 (|U of Mar|
|00002620| 79 6c 61 6e 64 2c 20 46 | 4c 42 29 22 3b 0d 0a 0d |yland, F|LB)";...|
|00002630| 0a 23 69 6e 63 6c 75 64 | 65 20 22 6d 61 74 2e 68 |.#includ|e "mat.h|
|00002640| 22 0d 0a 0d 0a 73 74 72 | 75 63 74 20 6d 61 74 72 |"....str|uct matr|
|00002650| 69 78 20 2a 0d 0a 6d 5f | 73 6f 6c 76 65 28 61 2c |ix *..m_|solve(a,|
|00002660| 20 62 29 0d 0a 72 65 67 | 69 73 74 65 72 20 73 74 | b)..reg|ister st|
|00002670| 72 75 63 74 20 6d 61 74 | 72 69 78 20 2a 61 2c 20 |ruct mat|rix *a, |
|00002680| 2a 62 3b 0d 0a 7b 0d 0a | 72 65 67 69 73 74 65 72 |*b;..{..|register|
|00002690| 20 73 74 72 75 63 74 20 | 6d 61 74 72 69 78 20 2a | struct |matrix *|
|000026a0| 61 5f 74 72 61 6e 73 2c | 20 2a 74 2c 20 2a 74 5f |a_trans,| *t, *t_|
|000026b0| 69 6e 76 2c 20 2a 74 32 | 2c 20 2a 78 3b 0d 0a 0d |inv, *t2|, *x;...|
|000026c0| 0a 69 66 20 28 28 61 2d | 3e 6d 5f 72 6f 77 73 29 |.if ((a-|>m_rows)|
|000026d0| 20 3c 20 28 61 2d 3e 6d | 5f 63 6f 6c 73 29 29 20 | < (a->m|_cols)) |
|000026e0| 7b 0d 0a 09 70 72 69 6e | 74 66 28 22 6d 5f 73 6f |{...prin|tf("m_so|
|000026f0| 6c 76 65 3a 20 74 6f 6f | 20 66 65 77 20 65 71 75 |lve: too| few equ|
|00002700| 61 74 69 6f 6e 73 5c 6e | 22 29 3b 0d 0a 09 72 65 |ations\n|");...re|
|00002710| 74 75 72 6e 28 4d 5f 4e | 55 4c 4c 29 3b 0d 0a 09 |turn(M_N|ULL);...|
|00002720| 7d 0d 0a 0d 0a 69 66 20 | 28 28 61 2d 3e 6d 5f 72 |}....if |((a->m_r|
|00002730| 6f 77 73 29 20 21 3d 20 | 28 62 2d 3e 6d 5f 72 6f |ows) != |(b->m_ro|
|00002740| 77 73 29 29 20 7b 0d 0a | 09 70 72 69 6e 74 66 28 |ws)) {..|.printf(|
|00002750| 22 6d 5f 73 6f 6c 76 65 | 3a 20 61 72 67 75 6d 65 |"m_solve|: argume|
|00002760| 6e 74 73 20 64 6f 6e 27 | 74 20 6d 61 74 63 68 3a |nts don'|t match:|
|00002770| 20 25 64 2c 20 25 64 2e | 5c 6e 22 2c 20 61 2d 3e | %d, %d.|\n", a->|
|00002780| 6d 5f 72 6f 77 73 2c 20 | 62 2d 3e 6d 5f 72 6f 77 |m_rows, |b->m_row|
|00002790| 73 29 3b 0d 0a 09 72 65 | 74 75 72 6e 28 4d 5f 4e |s);...re|turn(M_N|
|000027a0| 55 4c 4c 29 3b 0d 0a 09 | 7d 0d 0a 0d 0a 69 66 20 |ULL);...|}....if |
|000027b0| 28 62 2d 3e 6d 5f 63 6f | 6c 73 20 21 3d 20 31 29 |(b->m_co|ls != 1)|
|000027c0| 20 7b 0d 0a 09 70 72 69 | 6e 74 66 28 22 6d 5f 73 | {...pri|ntf("m_s|
|000027d0| 6f 6c 76 65 3a 20 61 72 | 67 32 20 6d 75 73 74 20 |olve: ar|g2 must |
|000027e0| 68 61 76 65 20 31 20 63 | 6f 6c 75 6d 6e 2e 5c 6e |have 1 c|olumn.\n|
|000027f0| 22 29 3b 0d 0a 09 72 65 | 74 75 72 6e 28 4d 5f 4e |");...re|turn(M_N|
|00002800| 55 4c 4c 29 3b 0d 0a 09 | 7d 0d 0a 0d 0a 61 5f 74 |ULL);...|}....a_t|
|00002810| 72 61 6e 73 20 3d 20 6d | 5f 74 72 61 6e 73 70 6f |rans = m|_transpo|
|00002820| 73 65 28 61 29 3b 09 09 | 2f 2a 20 41 27 20 2a 2f |se(a);..|/* A' */|
|00002830| 0d 0a 74 20 3d 20 6d 5f | 6d 75 6c 74 69 70 6c 79 |..t = m_|multiply|
|00002840| 28 61 5f 74 72 61 6e 73 | 2c 20 61 29 3b 09 09 2f |(a_trans|, a);../|
|00002850| 2a 20 41 27 20 41 20 2a | 2f 0d 0a 74 5f 69 6e 76 |* A' A *|/..t_inv|
|00002860| 20 3d 20 6d 5f 69 6e 76 | 65 72 74 28 74 29 3b 09 | = m_inv|ert(t);.|
|00002870| 09 09 2f 2a 20 28 41 27 | 20 41 29 2d 31 20 2a 2f |../* (A'| A)-1 */|
|00002880| 0d 0a 66 72 65 65 28 74 | 29 3b 0d 0a 69 66 20 28 |..free(t|);..if (|
|00002890| 74 5f 69 6e 76 20 3d 3d | 20 4d 5f 4e 55 4c 4c 29 |t_inv ==| M_NULL)|
|000028a0| 20 7b 0d 0a 09 70 72 69 | 6e 74 66 28 22 6d 5f 73 | {...pri|ntf("m_s|
|000028b0| 6f 6c 76 65 3a 20 6e 6f | 20 73 6f 6c 75 74 69 6f |olve: no| solutio|
|000028c0| 6e 5c 6e 22 29 3b 0d 0a | 09 72 65 74 75 72 6e 28 |n\n");..|.return(|
|000028d0| 4d 5f 4e 55 4c 4c 29 3b | 0d 0a 09 7d 0d 0a 74 32 |M_NULL);|...}..t2|
|000028e0| 20 3d 20 6d 5f 6d 75 6c | 74 69 70 6c 79 28 74 5f | = m_mul|tiply(t_|
|000028f0| 69 6e 76 2c 20 61 5f 74 | 72 61 6e 73 29 3b 09 2f |inv, a_t|rans);./|
|00002900| 2a 20 28 41 27 20 41 29 | 2d 31 20 41 27 20 2a 2f |* (A' A)|-1 A' */|
|00002910| 0d 0a 66 72 65 65 28 74 | 5f 69 6e 76 29 3b 0d 0a |..free(t|_inv);..|
|00002920| 66 72 65 65 28 61 5f 74 | 72 61 6e 73 29 3b 0d 0a |free(a_t|rans);..|
|00002930| 78 20 3d 20 6d 5f 6d 75 | 6c 74 69 70 6c 79 28 74 |x = m_mu|ltiply(t|
|00002940| 32 2c 20 62 29 3b 09 09 | 09 2f 2a 20 28 41 27 20 |2, b);..|./* (A' |
|00002950| 41 29 2d 31 20 41 27 20 | 42 20 2a 2f 0d 0a 66 72 |A)-1 A' |B */..fr|
|00002960| 65 65 28 74 32 29 3b 0d | 0a 0d 0a 72 65 74 75 72 |ee(t2);.|...retur|
|00002970| 6e 28 78 29 3b 0d 0a 7d | 0d 0a 2f 2f 67 6f 2e 73 |n(x);..}|..//go.s|
|00002980| 79 73 69 6e 20 64 64 20 | 2a 0d 0a 6d 61 64 65 3d |ysin dd |*..made=|
|00002990| 54 52 55 45 0d 0a 69 66 | 20 5b 20 24 6d 61 64 65 |TRUE..if| [ $made|
|000029a0| 20 3d 20 54 52 55 45 20 | 5d 3b 20 74 68 65 6e 0d | = TRUE |]; then.|
|000029b0| 0a 09 2f 62 69 6e 2f 63 | 68 6d 6f 64 20 36 34 34 |../bin/c|hmod 644|
|000029c0| 20 6d 5f 73 6f 6c 76 65 | 2e 63 0d 0a 09 2f 62 69 | m_solve|.c.../bi|
|000029d0| 6e 2f 65 63 68 6f 20 2d | 6e 20 27 09 27 3b 20 2f |n/echo -|n '.'; /|
|000029e0| 62 69 6e 2f 6c 73 20 2d | 6c 64 20 6d 5f 73 6f 6c |bin/ls -|ld m_sol|
|000029f0| 76 65 2e 63 0d 0a 66 69 | 0d 0a 2f 62 69 6e 2f 65 |ve.c..fi|../bin/e|
|00002a00| 63 68 6f 20 27 45 78 74 | 72 61 63 74 69 6e 67 20 |cho 'Ext|racting |
|00002a10| 6d 5f 74 72 61 6e 73 70 | 6f 73 65 2e 63 27 0d 0a |m_transp|ose.c'..|
|00002a20| 73 65 64 20 27 73 2f 5e | 58 2f 2f 27 20 3c 3c 27 |sed 's/^|X//' <<'|
|00002a30| 2f 2f 67 6f 2e 73 79 73 | 69 6e 20 64 64 20 2a 27 |//go.sys|in dd *'|
|00002a40| 20 3e 6d 5f 74 72 61 6e | 73 70 6f 73 65 2e 63 0d | >m_tran|spose.c.|
|00002a50| 0a 73 74 61 74 69 63 20 | 63 68 61 72 20 2a 73 63 |.static |char *sc|
|00002a60| 63 73 69 64 20 3d 20 22 | 40 28 23 29 6d 5f 74 72 |csid = "|@(#)m_tr|
|00002a70| 61 6e 73 70 6f 73 65 2e | 63 09 34 2f 35 2f 38 32 |anspose.|c.4/5/82|
|00002a80| 20 28 55 20 6f 66 20 4d | 61 72 79 6c 61 6e 64 2c | (U of M|aryland,|
|00002a90| 20 46 4c 42 29 22 3b 0d | 0a 0d 0a 23 69 6e 63 6c | FLB)";.|...#incl|
|00002aa0| 75 64 65 20 22 6d 61 74 | 2e 68 22 0d 0a 0d 0a 73 |ude "mat|.h"....s|
|00002ab0| 74 72 75 63 74 20 6d 61 | 74 72 69 78 20 2a 0d 0a |truct ma|trix *..|
|00002ac0| 6d 5f 74 72 61 6e 73 70 | 6f 73 65 28 6d 61 74 29 |m_transp|ose(mat)|
|00002ad0| 0d 0a 72 65 67 69 73 74 | 65 72 20 73 74 72 75 63 |..regist|er struc|
|00002ae0| 74 20 6d 61 74 72 69 78 | 20 2a 6d 61 74 3b 0d 0a |t matrix| *mat;..|
|00002af0| 7b 0d 0a 72 65 67 69 73 | 74 65 72 20 73 74 72 75 |{..regis|ter stru|
|00002b00| 63 74 20 6d 61 74 72 69 | 78 20 2a 72 65 73 75 6c |ct matri|x *resul|
|00002b10| 74 3b 0d 0a 72 65 67 69 | 73 74 65 72 20 69 6e 74 |t;..regi|ster int|
|00002b20| 20 72 6f 77 2c 20 63 6f | 6c 3b 0d 0a 0d 0a 6d 5f | row, co|l;....m_|
|00002b30| 63 72 65 61 74 65 28 72 | 65 73 75 6c 74 2c 20 6d |create(r|esult, m|
|00002b40| 61 74 2d 3e 6d 5f 63 6f | 6c 73 2c 20 6d 61 74 2d |at->m_co|ls, mat-|
|00002b50| 3e 6d 5f 72 6f 77 73 29 | 3b 0d 0a 0d 0a 66 6f 72 |>m_rows)|;....for|
|00002b60| 20 28 72 6f 77 20 3d 20 | 30 3b 20 72 6f 77 20 3c | (row = |0; row <|
|00002b70| 20 6d 61 74 2d 3e 6d 5f | 72 6f 77 73 3b 20 72 6f | mat->m_|rows; ro|
|00002b80| 77 2b 2b 29 0d 0a 09 66 | 6f 72 20 28 63 6f 6c 20 |w++)...f|or (col |
|00002b90| 3d 20 30 3b 20 63 6f 6c | 20 3c 20 6d 61 74 2d 3e |= 0; col| < mat->|
|00002ba0| 6d 5f 63 6f 6c 73 3b 20 | 63 6f 6c 2b 2b 29 0d 0a |m_cols; |col++)..|
|00002bb0| 09 09 6d 5f 76 28 72 65 | 73 75 6c 74 2c 20 63 6f |..m_v(re|sult, co|
|00002bc0| 6c 2c 20 72 6f 77 29 20 | 3d 20 6d 5f 76 28 6d 61 |l, row) |= m_v(ma|
|00002bd0| 74 2c 20 72 6f 77 2c 20 | 63 6f 6c 29 3b 0d 0a 0d |t, row, |col);...|
|00002be0| 0a 72 65 74 75 72 6e 28 | 72 65 73 75 6c 74 29 3b |.return(|result);|
|00002bf0| 0d 0a 7d 0d 0a 2f 2f 67 | 6f 2e 73 79 73 69 6e 20 |..}..//g|o.sysin |
|00002c00| 64 64 20 2a 0d 0a 6d 61 | 64 65 3d 54 52 55 45 0d |dd *..ma|de=TRUE.|
|00002c10| 0a 69 66 20 5b 20 24 6d | 61 64 65 20 3d 20 54 52 |.if [ $m|ade = TR|
|00002c20| 55 45 20 5d 3b 20 74 68 | 65 6e 0d 0a 09 2f 62 69 |UE ]; th|en.../bi|
|00002c30| 6e 2f 63 68 6d 6f 64 20 | 36 34 34 20 6d 5f 74 72 |n/chmod |644 m_tr|
|00002c40| 61 6e 73 70 6f 73 65 2e | 63 0d 0a 09 2f 62 69 6e |anspose.|c.../bin|
|00002c50| 2f 65 63 68 6f 20 2d 6e | 20 27 09 27 3b 20 2f 62 |/echo -n| '.'; /b|
|00002c60| 69 6e 2f 6c 73 20 2d 6c | 64 20 6d 5f 74 72 61 6e |in/ls -l|d m_tran|
|00002c70| 73 70 6f 73 65 2e 63 0d | 0a 66 69 0d 0a 0d 0a 0d |spose.c.|.fi.....|
|00002c80| 0a 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ca0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ce0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002cf0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d00| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002d90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002da0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002db0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002dc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002dd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002de0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002df0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e00| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002e90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ea0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002eb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ec0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ed0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ee0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ef0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f00| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f10| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f30| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f40| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002fa0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002fb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002fc0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002fd0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002fe0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ff0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+